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SECTION  1 

INTRODUCTION  AND  SUMMARY 

Hg-Cd-Te  is  a  semiconductor  material  particularly  useful  as 
an  infrared  detector.  Although  it  has  been  grown  for  some 
years,  there  are  still  problems  in  achieving  crystals  of  desired 
quality.  In  order  to  improve  the  crystal  growth  technique,  the 
knowledge  and  understanding  of  the  liquidus-solidus  phase  dia¬ 
gram  are  of  crucial  importance.  This  study  has  addressed  the 
theoretical  calculation  of  the  phase  diagram  of  HgCdTe. 

Heretofore,  we  have  known  that  the  pair  approximation  of 
the  cluster  variation  method  (CVM)  provides  a  good  theoretical 
phase  diagram  of  III-V  compound  semiconductors  wuch  as  GaAs  and 
In-As,  as  reported  in  Section  2  and  Appendix  A.  In  this  project 
we  have  applied  a  similar  theoretical  technique  to  Hg-Cd-Te, 
which  is  a  II-VI  semiconductor. 

The  special  qualitative  features  of  the  II-VI  liquidus  as 
compared  with  the  III-V  case  are  that  the  former  is  asymmetric 
with  respect  to  the  50-50  composition,  and  that  the  liquidus  is 
sharply  peaked  at  the  50%  composition.  In  order  to  explain 
these  two  features  theoretically  (Section  3  and  Appendix  B),  we 
have  allowed  vacancies  in  the  pseudo-lattice  model  of  the  liquid 
phase,  and  have  included  molecular  species  HgTe  and  CdTe  in  the 
liquid  phase.  It  has  been  discovered  that  the  peak  in  the 
liquidus  results  when  the  interaction  between  an  atomic  species 
(Hg,  Cd  or  Te)  and  a  molecular  species  (HgTe  or  CdTe)  is 
repulsive. 

Using  the  associated  liquid  model  developed  in  Section  3,  a 
theoretical  treatment  of  the  liquidus  of  the  ternary  Hg-Cd-Te 
has  been  developed  in  Section  4  (and  in  Appendix  C).  In  this 
treatment  no  lattice  defects  are  included  in  the  solid  phase. 

The  liquidus  in  both  the  Te-rich  region  and  that  in  the  Hg-rich 
region  agree  with  experiments  reasonably  well.  The  special 
feature  of  the  work  that  differs  from  the  previous  regular 


solution  theory  treatment  of  Brebnck  et  al.  is  that  the 
nearest-neighbor  interaction  parameters  used  in  formulating  the 
free  energy  of  the  liquid  phase  are  treated  as  being  independent 
of  temperature  and  composition. 

It  is  known  that  the  electrical  properties  of  the  semicon¬ 
ductor  depends  greatly  on  the  lattice  defects  in  the  Hg-Cd-Te 
crystal.  Thus,  in  Section  5  we  take  into  account  lattice  vacan¬ 
cies  and  antiatoms  in  the  crystal  (details  are  described  in 
Appendix  D) .  The  lattice  is  composed  of  two  fee  sublattices, 
one  for  Te  atoms  and  the  other  for  Hg  and  Cd.  The  theory  uses 
intersublattice  pair  probabilities  and  intrasublattice  pair 
probabilities  as  the  basic  variables  and  includes  intersublat¬ 
tice  and  intrasublattice  nearest-neighbor  interaction  energies. 

The  method  of  calculating  the  liquidus  and  solidus  using 
the  solid  with  defects  (of  Section  5)  is  explained  in  Section  6, 
and  the  details  are  presented  in  Appendix  B.  The  main  results 
are  the  number  of  Hg  vacancies  in  Figure  1(1  the  Hg  partial  pres¬ 
sure  in  Figure  11,  and  the  relation  between  the  two  in  Figure  12. 
These  results  agree  well  with  Vydyanathan ' s  experiments. 

Interstitial  atoms  have  not  been  taken  into  account  in 
Sections  5  and  6,  but  their  plausibility  is  discussed  in 
Section  7. 


. 


SECTION  2 


III-V  DIAGRAM  CALCULATIONS  USING  THE  "PAIR"  METHOD 

Calculations  of  phase  diagrams  of  III-V  systems  form  the 
basis  of  the  II-VI systems ,  which  are  the  main  concern  of  this 
study.  Therefore,  we  have  included  the  main  results  of  III-V 
calculations  in  Figures  1,  2,  and  3  even  though  they  had  been 
done  before  this  effort  was  begun.  They  were  calculated  by 
Stringfellow  and  Greene1  and  also  by  the  author.2  The 
formulation  in  Reference  2  is  basic  in  subsequent  work  developed 
in  this  study.  Therefore,  Reference  2  is  attached  as  Appendix  A 
in  this  report. 

In  modeling  and  calculating  the  free  energy  of  the  liquid 
phase,  a  pseudo-lattice  structure  is  assumed.  Configurations  of 
nearest-neighbor  pairs  on  the  lattice  are  chosen  as  the  basic 
variables,  and  thus  the  method  is  called  the  "pair”  method 
(which  is  also  called  the  quasi-chemical  approximation  or 
Be the ' s  approx imat ion ) . 

The  excellent  agreement  between  theory  and  experiments  in 
Figures  1,  2,  and  3  supports  the  usefulness  and  the  reliability 
of  the  pair  method.  Noteworthy  features  in  the  liquidus  curves 
in  Figure  1  are  that  the  curves  are  symmetric  with  respect  to 
the  50%  composition,  and  also  well  rounded  near  the  maxima. 
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Figure  1 .  Liquidus  curves  of 

III-V  semiconductors. 
Solid  curves  are 
theory  and  circles 
are  experiments. 


Figure  2.  The  calculated  iso¬ 
thermal  liquidus  of 
Ga-In-As  ternary  sys¬ 
tem  at  1 250  *K.  The 
numbers  along  the 
curve  are  x  in  the 
solid  composition, 
(GaAs) ,  (InAs)  . 
Four  tie  lines  are 
shown . 


Figure  3.  Pseudobinary  liquidus 
and  solidus  curves  of 
InAs-GaAs  (aj  and 
InSb-InAs  (b).  The 
solid  curves  are  cal¬ 
culations,  and  the 
circles  and  triangles 
are  experimental 
results. 


SECTION  3 

II-VI  BINARY  LIQUIDUS-SOLIDUS  CALCULATIONS  USING  THE 
"PAIR"  METHOD  AND  THE  ASSOCIATION  MODEL 

Different  from  the  III-V  cases  of  Figure  1,  the  liquidus 
curves  of  Hg-Te  and  Cd-Te  are  not  symmetric  at  right  and  left 
and  are  peaked  near  the  50%  composition.  In  order  to  take  into 
account  these  two  features/  the  model  of  the  liquid  phase  is 
modified  in  two  aspects:  vacancies  are  included  in  the  pseudo¬ 
lattice  structure  and  molecular  species  of  HgTe  or  CdTe  are 
added.  When  the  effective  interaction  of  molecular  species  and 
atomic  species  is  repulsive,  we  can  explain  the  sharp  peak  near 
the  50%  composition.  The  results  are  shown  in  Figures  4  and  5, 
and  are  published  in  Reference  3,  which  is  attached  in  this 
report  as  Appendix  B. 

The  model  of  the  liquid  state  in  which  the  molecular 
species  are  included  is  called  the  associated  model.  In  such  a 
model,  the  number  of  species  is  four  (i.e.,  Hg,  Te,  HgTe  and 
vacancy).  The  cluster  variation  method  (CVM)4  formulation  of 
the  pair  method  for  such  a  system  has  a  decisive  advantage  over 
the  quasi-chemical  approximation  treatment  of  the  same  problem, 
although  the  two  approaches  give  exactly  the  same  results, 
because  of  the  intuitive  ease  of  the  CVM  formulation. 
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Figure  4. 

Hg-Te  binary  liquidus. 
Dots  are  experiments. 
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SECTION  4 


Hg-Cd-Te  LIQUIDUS-SOLIDUS 


PHASE 


DIAGRAM 


In  the  binary  cases  of  Section  3,  the  main  input  to  the 
theory  are  the  values  of  interaction  potentials  between  nearest- 
neighbor  pairs  in  the  pseudo-lattice  structure  of  the  liquid 
phase.  These  same  interaction  potential  values  of  binary  cases 
are  used  for  the  ternary  liquid  phases  as  well. 

In  the  Hg-Cd-Te  liquid  phase,  we  work  with  three  atomic 
species,  two  molecular  species  and  vacancies  on  the  pseudo¬ 
lattice  structure.  The  details  are  in  Reference  5  which  is 
reproduced  as  Appendix  C.  The  main  results  are  presented  in 
Figures  6,  7  and  8,  which  indicate  good  agreement  between  theory 
and  experiments.  Figure  6  shows  the  Hg  corner  (of  the  triangu¬ 
lar  Gibbs  diagram),  and  Figure  7  and  8  depict  the  Te  corner.  As 
opposed  to  the  regular  solution  approach  which  uses  the  "point" 
variables  only  rather  than  the  "pair"  variables,  the  present 
method  treats  the  energy  parameters  as  being  independent  of 
temperature  and  composition,  and  common  to  both  the  Hg  and  Te 
corners. 
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FRACTION  OF  Cd 


Figure  7 


The  ternary  Hg-Cd-Te  liquidus  near 
the  Te  corner.  The  solid  curves 
are  the  theory,  and  the  broken  curves 
are  the  experiments.  The  numbers  0.2 
and  0.4  have  the  same  meaning  as 
Figure  6. 


FRACTION  OF  Cd  Te  IN  SOLID 


SECTION  5 


Hg-Cd-Te  SOLID  PHASE  WITH  LATTICE  DEFECTS 

In  the  treatments  of  Sections  3  and  4,  the  solid  phase  is 
regarded  as  perfectly  stoichiometric.  In  actuality,  the  mea¬ 
sured  hole  concentration  is  associated  with  Hg  vacancies.  It  is 
also  known  that  antiatoms  exist,  which  means  that  Te  atoms  may 
sit  on  the  Hg-Cd  sublattice,  and  Hg  or  Cd  atoms  may  sit  on  the 
Te  sublattice. 

The  structure  of  such  a  solid  phase  with  lattice  defects 
has  been  calculated  and  is  attached  as  Appendix  D.  The  Gibbs 
free  energy  for  the  solid  phase  in  Appendix  A  and  B  is  repre¬ 
sented  by  a  single  point  because  it  is  for  the  perfect  stoichi¬ 
ometric  state.  It  is  different  when  lattice  defects  are 
included.  Figure  9  shows  an  example  of  how  the  Gibbs  free 
energy  changes  as  a  function  of  the  vacancy  concentration  in  the 
Hg-Cd  sublattice. 


SECTION  6 


Hg-Cd-Te  LIQUIDUS  AND  SOLIDUS  INCLUDING 
LATTICE  DEFECTS  IN  SOLID 


Combining  the  solid  phase  theory  in  Section  5  (and  Appendix 
D)  and  the  liquid  phase  theory  used  previously  in  References  3 
and  5,  we  have  formulated  the  coexistence  of  the  liquid  and  solid 
phases.  Coexistence  is  derived  by  using  the  conditions  in  which 
the  chemical  potentials  of  each  species  in  the  liquid  and  solid 
phases  are  equal  and  when  the  two  phases  are  under  the  same 
pressure. 

The  details  of  the  formulation  are  described  in  Appendix  E. 
The  main  results  are  shown  in  Figures  10,  11  and  12.  In  calcu¬ 
lating  these  figures,  the  interatomic  potentials  are  adjusted  so 
that  we  can  come  to  reasonable  agreement  with  experimental  data 
by  Vydyanath6 .  These  curves  are  calculated  for  the 
Hg0.8Cd0.2Te  solid  composition. 

Figure  10  plots  twice  the  number  of  vacancies  in  the  Hg-Cd 
sublattice  per  cm3  against  1/T  for  the  liquid-solid  coexistence 
condition.  The  solid  curves  are  the  present  theory  and  the 
broken  curves  are  the  experimental  hole  concentration  due  to 
Vydyanath6.  The  upper  curves  are  when  the  solid  coexists  with 
the  Te-rich  liquid,  and  the  lower  curves  are  for  the  Hg-rich 
liquid.  The  fair  agreement  between  the  theory  and  the  experi¬ 
ments  supports  the  view  that  a  Hg  vacancy  contributes  two  holes 
in  the  valence  bond6. 

In  adjusting  energy  parameters  to  change  the  general  shape 
of  the  curves,  we  note  two  properties.  Notations  are  from 
Appendix  E. 

[I]  As  we  make  ellf  e22.  and  e12  more  negative,  keeping 

el2  5  £ 1 2  “  Ui  l  +e22  fl-xed'  the  2  vHg  curves  in  Figure 
10 come  down  parallel  to  themselves. 

[II]  As  we  decrease  ult  and  u22  together,  the  vHg  curve 
for  the  Te-rich  case  does  not  change,  while  the  vHg  for 

the  Hg-rich  case  increases. 
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Figure  11 


.  Hg  partial  pressure  plotted  against 
Solid  curves  are  the  present  theory 
the  broken  curves  are  Vydyanath's 
experiments. 6 


Hg  PARTIAL  PRESSURE  PHg  ATM 

The  theoretical  relation  between  Hg  vacancies 
and  the  partial  pressure  of  Hg.  The  two  out¬ 
side  curves  are  for  liquid-solid  coexistence 
conditions,  and  the  line  connecting  them  is 
for  states  of  intermediate  Hg  pressure. 


The  partial  pressure,  PHg,  of  Hg  in  the  vapor  phase 
can  be  calculated  from  the  chemical  potential  of  Hg  when  we 
assume  that  Hg  in  the  vapor  phase  behaves  as  an  ideal  gas. 

Figure  11  plots  PHg  against  1/T.  The  solid  curves  are  the 
present  theory  and  the  broken  curves  are  Vydyanath ' s 
experiments6 .  The  theory  uses  the  same  energy  parameters  as 
those  for  Figure  10.  The  PHg  value  at  500 #C  for  the  Hg- 
rich  growth  is  chosen  at  the  experimental  value,  6.5  atom. 

In  adjusting  the  theoretical  curves,  we  noted  the  following 
property: 

[III]  The  energy  parameter,  e33,  controls  the  direction  of 
the  two  PHg  curves.  When  £33  is  more  negative,  the  PHg  curves 
point  further  upward. 

The  energy  parameter  values  that  we  used  are  in  Table  1 . 

The  energy,  u^j,  is  for  intersublattice  pairs.  These  pairs  are 
predominantly  Hg-Te  and  Cd-Te,  and  thus  the  u^j  values  for  these 
pairs  are  negative  (attractive)  while  the  rest  are  positive 
(repulsive).  The  is  for  intrasublattice  pairs?  its  sign  is 
opposite  to  that  of  u^j.  The  values  for  the  Hg-Te  and  Cd-Te 
pairs  were  chosen  as  zero,  but  they  can  be  positive;  since  the 
number  of  these  pairs  are  very  small,  the  actual  values  are  not 
important  in  the  phase  diagram  calculations. 

After  a  crystal  is  made  from  a  Hg-rich  melt  or  a  Te-rich 
melt,  the  lattice  defect  properties  can  be  changed  by 
controlling  the  Hg  partial  pressure  on  the  crystal.  Figure  12 
shows  the  relation  between  2v..  and  P„  .  These  intermediate 
points  are  calculated  by  changing  the  chemical  potential  value 
in  the  theory.  The  intermediate  points  lie  almost  on  a  line  in 
Figure  12,  in  a  manner  similar  to  the  experiments  reported  by 
Vydyanath  6 . 


SECTION  7 


INTERSTITIAL  Hg  IN  THE  Hg-Cd-Te  CRYSTAL 

In  the  treatment  of  Sections  5  and  6,  we  neglected  the 
possibility  of  interstitial  atoms.  Since  Hg  is  volatile,  it  is 
worthwhile  to  consider  interstitial  Hg  atoms. 

As  is  discussed  in  the  Introduction  (and  Figure  1 )  of 
Appendix  D,  the  Hg-Cd-Te  crystal  is  made  of  two  fee  sublattices. 
We  call  the  white  fee  sublattice  in  Figure  1  of  Appendix  D  as 
I,  and  the  black  sublattice,  II.  When  we  shift  the  fee  I  and 
fee  II  by  a  half  of  the  cube  edge,  they  occupy  interstitial 
positions.  We  can  also  call  the  new  fee  sublattices  III  and  IV. 
The  small  corner  cube  shown  by  dotted  lines  in  Figure  1  of 
Appendix  D  now  has  the  structure  shown  in  (a)  of  Figure  13,  and 
an  adjacent  small  cube  becomes  (b).  There  are  also  small  cubes 
of  structures  in  Figure  13  (c)  and  (d). 

Finally,  we  can  consider  whether  a  Hg  atom  can  easily 
assume  an  interstitial  position  II  or  IV.  In  Figure  13,  we  see 
that  a  I  site  and  a  III  site  have  the  same  nearest  neighbor 
configuration,  while  a  II  and  a  IV  site  have  nearest  neighbors 
of  their  own.  Since  in  the  original  structure,  Hg  atoms  sit 
predominantly  on  I  sites,  we  see  that  it  is  quite  easy  for  a  Hg 
atom  to  come  to  a  III  interstitial  site.  It  is  therefore  rea¬ 
sonable  to  take  into  account  in  the  theory  the  interstitial  Hg 
atoms  sitting  on  the  III  interstitial  sublattice  sites.  The 
number  of  such  interstitial  Hg  atoms  will  depend  on  the  Hg-Hg 
and  Hg-Cd  interactions  between  I-III  sites.  Note  that  the  I-III 
distance  is  larger  than  the  I-II  distance,  but  smaller  than  the 
intrasublattice  I-I  distance. 

Such  intersitial  atoms  have  not  been  treated  in  the  present 
report,  but  are  to  be  included  in  future  studies. 
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Liquidua-iolidul  phase  diagrams  of  m-V  semiconductors  are  calculated  using  the  Cluster  Variation  Method  and  related 
techniques  developed  recently.  The  liquid  phase  is  treated  assuming  a  lattice  model.  An  equilibrium  state  is  derived  by 
minimizing  the  grand  potential  (d),  keeping  the  chemical  potential  fixed.  The  resulting  nonlinear  equations  are  solved 
using  the  Natural  Iteration  Method  which  has  a  property  that  the  value  of  6  always  decreases  at  each  iteration  cycle. 

The  phase  boundary  is  derived  from  an  intersection  of  two  <S  curves;  the  correct  forms  of  <j  values  to  be  used  are 
derived  bom  geometrical  considerations.  Ternary  phase  diagrams  calculated  for  In-Ga-As  and  In-Sb-As  agree  well  with 
previous  calculations  of  StringfeOow  and  Greene  and  with  experimental  data. 

For  ternary  cases,  a  tie  line  connecting  coexisting  liquid  and  solid  phases  are  proved  to  be  orthogonal  to  the  phase 
boundary  curve  in  the  chemical  potential  diagram. 


All  existing  theoretical  treatments  of  phase 
diagrams  (liquidus  and  solidus)  of  I1I-V  and 
II-VI  compounds  make  use  of  Vieland’s  work 
[5],  which  establishes  the  difference  between 
reference  levels  of  the  liquid  and  the  solid  phases 
of  the  same  material.  We  will  also  use  the  rela¬ 
tion  developed  by  Vieland. 

In  formulating  the  liquidus-solidus  diagram 
the  free  energy  of  the  solid  is  relatively  easy  to 
figure  out.  However,  to  treat  the  liquid  phase 
accurately  is  a  difficult  task  and  one  beyond  the 
scope  of  the  present  work.  Therefore,  we  follow 
S-G  and  use  a  lattice  model  for  the  liquid  phase. 
Since  the  formulation  is  done  using  the  pair 
approximation  (which  is  equivalent  to  the  QCM 
used  in  ref.  1)  of  the  CVM,  we  do  not  need  to 
specify  the  lattice  structure  for  the  liquid  phase 
except  the  coordination  number  z. 

2.  Sketch  of  S-G’s  work  of  binary  liquids 

This  section  reviews  those  formulations  of  the 
S-G  paper  [1]  that  are  relevant  to  our  develop¬ 
ment  in  later  sections  in  order  to  show  the  rela¬ 
tions  between  their  work  and  the  present  one. 
We  do  not  use  their  notation,  but  follow  the 
standard  formulation  of  the  pair  approximation 
of  the  CVM. 

In  treating  the  liquid  phase,  S-G  use  the  lat- 


Some  years  ago  Stringfellow  and  Greene  [1], 
to  be  cited  as  S-G,  applied  the  pair  ap¬ 
proximation  of  the  quasichemical  method 
(QCM)  (2]  to  calculate  liquidus  and  solidus  phase 
diagrams  of  m-V  binary  and  ternary  semicon¬ 
ductor  systems.  Their  results  for  In-Ga-As  and 
In-As-Sb  agree  very  well  with  experiments. 

After  their  work,  the  present  author 
developed  a  new  method  of  calculating  phase 
diagrams  [3,4]  as  an  application  of  the  hierar¬ 
chical  cluster  variation  method  (CVM).  The  new 
method  in  refs.  3  and  4  uses  the  grand  potential 
minimization  rather  than  matching  of  individual 
chemical  potentials  as  has  been  done  by  S-G 
and  others,  and  also  uses  the  natural  iteration 
method  (NIM)  in  solving  simultaneous  nonlinear 
algebraic  equations.  The  new  approach  is  more 
compact  than  other  existing  methods  and  thus 
presents  an  easier  means  of  generalizing  it  to 
more  complicated  systems. 

The  present  paper  reports  how  the  new 
method  is  formulated,  using  the  systems  treated 
by  S-G  as  examples.  The  numerical  results  are 
exactly  the  same  as  those  of  S-G,  but  some  of 
the  formulae  are  interpreted  differently  and  also 
additional  properties  of  the  systems  not  reported 
by  them  are  presented. 
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tice  model  of  die  coordination  number  z  *  2*> 
(since  we  often  use  zt 2,  it  is  written  as  u  in  this 
paper).  Each  lattice  point  is  occupied  by  either 
an  A  atom  (designated  by  a  subscript  i m  1)  or  a 
B  atom  (i  *  2),  and  we  do  not  allow  vacancies  in 
the  lattice.  The  probability  of  finding  an  <th 
species  on  a  lattice  point  is  written  as  x,  (i  *  1  or 
2).  The  probability  of  finding  an  i-j  nearest- 
neighbour  pair  is  written  as  y#  (i,  /  —  1  or  2). 
These  variables  satisfy  the  mutual  relations 


land  2, 


and  the  symmetry  relation 

y«  *  yu.  (22) 

The  normalization  relation  is 

x.  +  zj-l.  (2.3) 

The  potential  energy  of  the  lattice  liquid  is 
written  as 

E«2«/V  (2.4a) 

where  2N  is  the  total  number  of  lattice  points  in 
the  system,  and  is  the  potential  energy  for  a 
nearest-neighbor  pair  i-j.  The  energy  expression 
(2.4a)  can  be  rewritten  using  (2.1)  and  (2.2)  as 

E  -  2atN («„xi  +  «axi  -  4eyu),  (2.4b) 

where  c  is  defined  as 

4*  —  «u  +  —  2«i2-  (2.5) 

The  entropy  according  to  the  pair  ap¬ 
proximation  of  QCM  is  equivalent  to  that  of 
CVM  and  is  written  as 

S  -  2 kN  £(2 at  - 1)  x,  In  x,  -  to  yv  In  yv  j, 

(2.6) 

where  k  is  the  Boltzmann  constant.  Using  E  in 
(2.4b)  and  S  in  (2.6),  we  can  write  the  free 
energy 

F-E-TS  (2.7) 

as  a  function  of  x  and  y  values.  In  (2.7),  T  is  the 
absolute  temperature.  Note  that  5  in  (2.6)  is  for 
2 N  lattice  points. 


When  the  x  values  are  fixed,  there  is  only  one 
independent  variable  for  which  we  can  choose 
yu.  The  equilibrium  state  is  derived  by  minimiz¬ 
ing  F  in  (2.7)  with  respect  to  yu  keeping  T  and 
the  x  values  fixed.  The  differentiation  leads  to 

yuya-e-^yk,  (2.8) 

where  we  define 

0  - 1  IkT.  (2.9) 

Eq.  (2.8)  corresponds  to,  for  example,  (3a)  of 
S-G,  and  we  can  derive  the  identity  between 
S-G’s  /Jab  and  our  e  as 


Eq.  (2.8)  can  be  solved  with  the  aid  of  (2.1)  as 

yiz  *  ya  “  2xiXj/(1  +  r),  (2.11) 

where  r  is  written  as  0  in  S-G  and  is  defined  as 

r  -  (1  +  4x,x2(e-4*  - 1)1"*.  (2. 12) 

In  deriving  the  chemical  potentials,  it  is  helpful 
to  use 

Ni»2NXh  i -land  2,  n 

N¥-2 JVy*  4/ -land  2.  ^  ; 

In  terms  of  N  values,  we  can  write  F  in  (2.7) 
explicitly  as 

F  *  «(«nN|  +  CaNj  —  4 eNu) 

-kT  [ (2»  - 1)  T  2(N,)  -  m  2  Z  -W) 

-(«-l)if(2JV)],  (2.14) 

where  the  if  function  is  defined  as 

j?(x)"xlnx-x.  (2.15) 

When  F  is  a  minimum  with  respect  to  y^  = 
S\jj(2N)  and  (2.8)  holds,  the  chemical  potentials 
ft,  are  defined  from  (2.14)  by  differentiation  as 

*  ■  (m)v„  -  ■ + 1 W  [>n *y +  « in  ] 
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S-G  calls  the  last  ratio  the  activity  coefficients: 

*•[$  r and  »-&r  (2i?) 

Using  the  yu  expression  in  (2.11),  we  can  verify 
that  (2.17)  agrees  with  S-G’s  eq.  (4). 

Note  that,  in  eqs.  (2.8)  and  (2.12),  the  only 
energy  parameter  we  need  for  the  equilibrium 
properties  of  the  system  is  «,  defined  in  (2.5).  The 
individual  eu  and  «B  appear  in  (2.16),  but  they 
play  the  role  of  only  defining  the  zero  levels  of 
the  chemical  potentials 

3.  Grand  potential  formulation  of  the  liquid 
phase 

In  S-G’s  formulation  (review  above)  and  in 
most  current  formulations  of  phase  equilibria, 
one  starts  with  the  free  energy  and  minimizes  it 
keeping  the  composition  fixed  in  order  to  obtain 
the  equilibrium  state.  The  grand  potential  for¬ 
mulation  to  be  presented  in  this  section  leads  to 
exactly  the  same  results  as  those  of  the  previous 
section,  but  has  several  advantages  in  mathema¬ 
tical  convenience. 

We  start  with  the  grand  potential,  which  is 
defined  as 

(3  -F-2AM,  (3.1) 

I 

where  fi,  are  related  to  the  chemical  potential 
but  for  a  moment  can  be  regarded  as  Lagrange 
multipliers  for  adjusting  the  number  of  atoms  S,. 

To  facilitate  our  further  discussions,  we  write 
<j  in  (3.1)  in  full  using  (2.4a)  and  (2.6): 

d>-0(3/(2N) 

■  - 1)  [S  &(*<) + 2  ■£(*>)] 

+“X%)~("'i) 

-  &  [£  xA  +  £  +0A  (l  -  £  y,),  (3.2) 

where  fik  is  a  Lagrange  multiplier  used  for  the 
normalization  of  yv: 

1  -£*/•  (3.3) 


The  if  function  is  defined  in  (2.15).  In  (3.2),  the 
Ziffo)  sum  and  sum  are  repeated  and 

halved  for  the  sake  of  later  mathematical  con¬ 
venience. 

To  derive  the  equilibrium  state,  we  minimize 
the  grand  potential  <5  in  (3.2)  with  respect  to  the 
pair  variables  y^  treating  all  four  of  them  as 
being  independent  (in  contrast  to  treating  only 
yi2  as  independent  in  section  2)  and  keeping  /2,  as 
fixed  (rather  than  keeping  x,  as  fixed,  as  in 
section  2).  In  so  doing,  we  regard  x,  and  x,  in 
(3.2)  as  different  linear  combination  of  y„: 

Xi  «  2  y«  and  xi  *  2  y*  (3.4) 

/  * 

When  we  minimize  <P  in  (3.2),  keeping  0  and  /l, 
fixed,  we  arrive  at 

-  H2»  - 1)  Into; ) + ft  In  y„ 

-i0(/L  +  6)-A0*O  (3.5) 

Or 

yv  *  ea,*'to;)°”*,w2“) 

xexp[— +  #*»)]•  (3  6) 

We  now  see  that  the  normalization  (3.3),  the 
reduction  relations  (3.4),  and  the  “superposition" 
relations  (3.6)  form  a  set  of  equations  to  be 
solved  for  y*’s.  The  N1M  solves  them  systema¬ 
tically  as  follows. 

Note  that,  for  the  NIM.  we  fix  the  /3  and  £'s. 
We  start  with  the  first  input  values  x,  and  x:. 
Different  from  the  Newton-Raphson  method 
(NRM),  any  guess  values  for  the  first  input  lead 
to  a  converged  solution.  The  initial  guess  value  is 
used  on  the  right-hand  side  of  (3.6)  without  the 
normalization  factor  exp(A $/<»),  to  derive  the 
first  output  y#  exp(— A (31  at).  These  quantities  are 
used  in  (3.3)  to  obtain  the  normalization  factor 
exp(A0/w),  and  hence  the  output  y„  in  (3.6).  This 
set  of  the  first  output  yv  is  used  on  the  right-hand 
side  of  (3.4)  to  derive  the  second  input  set  x,  and 
Xj.  The  iteration  cycle  is  repeated  until  the  con¬ 
secutive  outputs  are  sufficiently  dose  to  each 
other. 

It  was  proved  previously  [6]  that  the  grand 
potential  <t>  in  (3.2)  atways  decreases  step  by  step 
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at  each  iteration  cycle.  The  proof  is  reviewed  in 
the  appendix.  This  proof  guarantees  that  the 
iteration  always  converges  to  a  minimum  of 
When  $  is  a  minimum  and  the  set  of  equations 
(3.3X  (3.4)  and  (3.5)  is  satisfied,  we  can  simplify 


the  expression  of  ♦  by  forming 

(3.7a) 

or 

A-  — 

A  IN' 

(3.7b) 

Thus,  the  Lagrange  multiplier  introduced  in  (3.2) 
for  the  normalization  of  yv’s  has  the  meaning  of 
die  grand  potential  per  lattice  point  in  the  equil¬ 
ibrium  state. 

One  particular  case  of  present  interest  to  us  is 
when  i  -  j  in  (3.5),  which  leads  to 

A  +  A  «■  +  kT  £ln  x,  +  at  In  (3.8) 

Compared  with  (2.16)  we  see  that  A  of  the 
present  section  is  related  to  the  chemical  poten¬ 
tial  Mi  derived  in  section  2  as 

M<”M<  +  A.  (3.9) 

Because  A  is  only  an  additional  constant,  we  can 
safely  also  call  A>  used  in  (3.1)  and  (3.2),  a 
chemical  potential. 

Use  of  the  meaning  of  A  in  (3.7b)  and  the 
relation  for  A  in  (3.9)  leads  to  a  useful  relation 
from  the  6  expression  (3.1): 

F-£m*Ni.  (3.10) 

Note  that  in  this  treatment  of  the  liquid  phase 
there  is  no  difference  between  the  Helmholtz 
free  energy  and  the  Gibbs  free  energy  because 
the  volume  is  kept  fixed  in  the  model. 

In  the  previous  section  it  was  shown  that  when 
the  composition  xt  (and  x:)  and  the  temperature 
T  are  given,  die  equilibrium  values  of  yu,  ytt, 
and  yn  can  be  determined  by  the  parameter  « 
only,  as  we  see  in  (2.11)  and  (2.12).  Therefore, 
we  expect  that  the  yv's  in  (3.6)  should  also 
depend  only  on  c  and  not  «tt  and  «a  explicitly, 
although  the  latter  appear  in  (3.6).  We  can  derive 


the  formulae  for  y^’s  which  depend  only  on  «  as 
follows. 

We  write 


y«-*-i£(y¥  +  y#)  (3.11) 

and  use  this  in  the  energy  expression  in  (2.4a): 

E  *  **  {&  *  [*  -  *5  <* + » >] + & l  *■*  1 

+«#)]>«]. 

(3.12) 


We  define 


+ <#). 

which  can  be  written  explicity  as 

•u  "  in  m  0, 

“  *2i  “  “2e, 


(3.13a) 

(3.13b) 


where  c  is  defined  in  (2.5).  When  (3.13a)  is  used 
in  (3.12X  we  can  write 

E-2wN  +  ^  (3- 14) 

When  this  expression  of  energy  is  used  in  (3.2), 
two  changes  need  to  be  made:  e«  should  be 
replaced  by  defined  in  (3.13),  and  A  should 
be  replaced  by 

AmA~ i  m  1  and  2.  (3.15) 


These  two  changes  appear  in  the  minimization 
relations  in  (3.5)  and  (3.6). 

In  the  NIM,  to  solve  the  y/s  from  (3.6),  which 
is  now  modified,  together  with  the  subsidiary 
conditions  (3.3)  and  (3.4),  we  fix  the  value  of  A< 
which  now  replaces  A  in  (3.6),  and  do  the  itera¬ 
tion.  The  A’s  are  the  parameters  that  control  the 
concentration  x*  and  the  fact  that  A  is  made  up 
of  two  terms  A  and  in  (3.15)  is  immaterial  as 
far  as  the  equilibrium  solution  is  concerned.  It  is 
to  be  noted  that,  since  <v  in  (3.13)  depends  on  e 
only  and  not  on  <n  and  *n  separately,  the  equil¬ 
ibrium  solution  for  the  y/s  depends  on  <  only  in 
agreement  with  what  we  noted  at  the  end  of 
section  2.  The  relations  between  the  A  and  x,’s 
and  A  are  also  independent  of  eu  and  By 
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using  the  modified  equation  (3.5),  modified  so 
that  <0,  ft*  and  it,  are  replaced  by  «#,  £,  and  /2„ 
and  by  comparing  with  the  relations  for  the  jt/s 
in  (2.16),  we  can  derive 

ft,  ■  fi,  +  A  +  a**,  i  ■  1  and  2,  (3-16) 

which  agrees  with  (3.9)  when  (3.15)  is  used.  The 
relations  in  (3.16)  agree  with  the  fact  that  «u  and 
ea  appear  explicitly  when  we  write  the  ft,  values 
in  (2.16). 

The  actual  numerical  computation  of  the 
equilibrium  state  is  made  with  iv  defined  in 
(3.13)  replacing  <#,  and  with  ft,  replacing  it,  in 
(3.6).  Since  we  assume  that  there  are  no  vacan¬ 
cies  in  the  system,  we  have 

X1  +  X2-I.  (3.16) 

This  imposes  a  mutual  relation  between  ft,  and 
(t 2,  and  thus  there  is  one  degree  of  freedom  in 
choosing  them.  In  calculating  binary  systems  we 
usually  treat  the  two  components  symmetrically 
and  use  the  condition 

fit  +  fi- 2*0.  (3.17) 

Sometimes,  however,  we  use  the  condition 
/*2-0,  (3.18) 

as  we  see  in  section  9. 

When  we  combine  (2.10)  with  (3.13b)  we  see 
that  our  in  is  related  to  of  S-G  [1]  as 

in  m  <21  *  ■OW*.  (3.19) 

where  z  is  the  coordination  number  used  for  the 
liquid  phase. 

4.  Derivation  of  the  binary  liquid  us 

The  liquidus  of  a  ITI-V  system  like  Ga-As, 
In- As,  Ga-Sb,  and  In-Sb  schematically  looks 
like  fig.  la.  The  point  P  is  in  the  liquid  phase, 
and  (?  is  on  the  solid  phase  of  the  composition 
AB.  The  width  of  the  solid  phase  can  be  neglec¬ 
ted  in  this  study  and  thus  we  assume  that  the 
solid  phase  has  the  fixed  composition  X|  -  x2  ■  i 
Since  the  two  points  P  and  Q  coexist,  the 
tangent  of  the  F  curve  at  P  in  fig.  lb  passes 
through  Q,  which  is  the  value  of  the  free  energy 
in  the  solid  phase.  This  condition  is  now  expres- 


Fi|.  1.  Schematics  of  the  liquidus  curve  (a)  and  the  free 
energy  F  (b)  for  a  binary  semiconductor.  Points  P  and  Q 
coexist. 

sed  in  mathematical  form.  First,  the  equation  of 
a  tangent  of  the  F  curve  at  P  can  be  written  as 

F_^(-)W).  (4.1) 

Using  the  normalization  of  x,  in  (2.3)  and  the 

definition  of  ft,  in  (2.16),  we  can  write 

G£r-[(©.-(©X 

-[ftp- finis.  (4.2) 

When  (4.2)  and  (3.10)  are  substituted  into  (4.1), 
we  obtain 

F~[fi\”Xl  +  tiPxi]2N.  (4.3) 

This  is  the  equation  for  the  tangent  to  the  F 

curve  in  fig.  lb  going  through  the  point  P  Eq. 
(4.3)  looks  similar  to  the  genera!  expression  of  F 
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in  (3.10),  but  is  different  in  the  sense  that  xx  and 
x2  are  variables  independent  of  the  point  P. 

The  next  important  step  is  to  estimate  the 
value  of  the  free  energy  of  the  solid  phase  at  the 
point  Q.  For  this  purpose  we  follow  S-G  [1]  and 
use  the  relation  originally  derived  by  Vieland  [5]. 
Vieland’s  relation  says  that  the  free  energy  of 
one  mole  of  molecular  AB  in  the  solid  phase  is 
written  as 

Fi&  -  F™  ~  (Tm  -  T)ASm,  (4.4) 

where  Fi<r>  is  the  free  energy  of  the  supercooled 
liquid  system  at  Q'  in  fig.  lb.  This  state  is 
“supercooled”  because  it  is  not  stable  compared 
with  the  solid  phase.  In  (4.4)  Tm  is  the  melting 
temperature  of  the  solid  AB,  as  indicated  in  fig. 

la,  and  AS.  is  the  entropy  of  melting  of  one 
mole  of  the  solid.  Eq.  (4.4)  is  important  in  fixing 
the  relative  zero  points  of  the  solid  and  liquid 
phases. 

We  can  use  (3.10)  in  writing  F0'’, 

(4.5) 

where  2N  is  the  total  number  of  atoms  in  the 
liquid  phase,  and  N,  =»  -  N. 

In  the  (F,  x2)  plane,  a  line  passing  through  the 
solid  point  (F0),  x<0))  is  written  as 

F-F«» -«(*!- **<”),  (4.6) 

where  a  can  be  any  value  and,  as  shown  in  fig. 

lb, 

xf'-l  (4.7) 

Eq.  (4.6)  is  to  be  compared  with  the  liquid  coun¬ 
terpart  (4.1). 

The  condition  that  the  solid  at  Q  in  fig.  lb  is  in 
equilibrium  with  the  liquid  at  P  is  that  the 
tangent  lines  (4.1)  and  (4.6)  become  identical. 

To  put  this  in  mathematical  form,  it  is  con¬ 
venient  to  introduce  the  following  two  expres¬ 
sions: 

«.8, 

<5,-FJ°>-oxi°».  (4.9) 

Since  no  confusion  is  expected  with  the  analysis 
in  section  3,  we  may  call  (5,  and  6,  grand  poten¬ 


tials  also.  In  (4.9),  x[Q)  actually  is  i,  as  we  see  in 
(4.7).  Since  a  in  (4.6)  and  hence  in  (4.9)  can  be 
any  value,  and  since  we  require  that  the  lines  in 
(4.1)  and  (4.6)  coincide,  we  require 


a 


fd F\n 
\dxj  ' 


(4.10) 


When 

<5.  -  <J„  (4.11) 

we  can  prove  that  the  tangent  to  F  at  P  goes 
through  the  solid  point  Q  as  follows.  Using  G,  in 
(4.8),  we  can  write  eq.  (4.1)  for  the  tangent  to  the 
F  curve  at  P  as 


We  can  show  using  (4.10)  that  the  point 
(Fw*,xfn)  lies  on  this  line  when  (4.11)  holds, 
because 

(4.13) 

For  numerical  computation,  (4.11)  can  be  fur¬ 
ther  simplified.  We  use  relations  (4.2),  (4.3), 
(4.4),  and  (4.5)  to  obtain 

<5,-2 NpP  (4.14a) 

<5,  -  AfOsP  +  M°V (Ta  -  T)ASm 

(4.14b) 

Substituting  thes  two  expressions  into  (4.11) 
yields 

N0*r  +#*?’)  -  N(mP  +  + (Tm  -  T)ASm. 

(4.15) 

In  reducing  the  equation  in  the  next  step,  we 
have  to  be  careful  about  the  mole  number.  When 
one  mole  of  GaAs  melts,  it  produces  one  mole  of 
Ga  and  one  mole  of  As.  Thus,  the  total  number 
of  atoms  in  the  liquid  phase  is  two  moles.  If  ASm 
in  (4.15)  is  measured  for  one  mole  of  GaAs.  S  in 
(4.15)  is  the  Avogadro  number  NAwf.  Thus,  by 
dividing  both  sides  by  NAvt  we  obtain 

Mr  +  nr  -  MP  +  MP -  (Tm  -  T)  (4. 16) 

4^Avg 


We  can  further  simplify  this  by  using  (3.16)  and 
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Table  I 

Parameter!  wed  in  the  calculation  of  Ga-As,  In-Aa  and 
In-Sb  binary  liquidus. 


Tm  (K)  AS.  i„  (cal/mok) 

(cal  deg'1  mole'? 


GaAs 

1511 

14.7 

-4380/6 

InAs 

1215 

14.7 

-6070/6 

InSb 

803 

13.3 

-3980/6 

*In  table  1  of  S-G  [1]  the  unit  of  AS.  is  written  sseu.  Usually 
eu  means  entropy  unit  and  represents  the  number  in  units  of  the 
pa  constant  If.  Their  eu  should  read  cal  deg*1  mote"1.  T.  is  the 
melting  temperature  of  the  compound.  AS.  is  the  entropy  of 
melting  of  the  compound,  and  « is  the  energy  parameter  of  the 
pair  in  the  liquid  phase  in  eq.  (2.5s). 

adopting  the  condition  (3.17)  as 
A(i',,A(07_(T.-T)^-.  (4.17) 

The  numerical  calculation  of  finding  the  liqui- 
dus  at  a  temperature  T  is  done  in  the  following 
way.  We  use  the  values  of  ASm  and  Tm  given  in 
table  I,  which  shows  the  values  used  by  S-G  [1], 
First,  we  solve  the  supercooled  liquid  state  for 
the  temperature  and  determine  the  value  A (<n  for 
the  point  O'  in  fig.  lb;  for  that  purpose,  we  use 
the  NIM  of  section  3. 

Next  we  vary  faf*  in  the  liquid  phase  and  find  the 
value  A  ^  which  satisfies  (4.17).  Then  the  point  P  is 
on  the  liqukhis.  We  have  used  this  formulation  and 
calculated  the  liquidus  curves  for  GaAs,  InAs,  and 
InSb.  The  results  agree  exactly  with  those  of  fig.  3 
in  S-G  [1]  and  thus  fit  nicely  with  experiments.  To 
make  the  presentation  complete,  we  plot  the 
curves  in  fig.  2. 

The  C|2  values  in  table  I  are  negative.  This 
shows  from  the  definition  of  ««  in  (3.13)  that  the 
different  species  of  the  pair  in  the  table  have  the 
tendency  to  attract  in  the  liquid  phase. 

5.  Ternary  liquid  phase 

In  section  3  we  worked  out  the  binary  liquid 
phase  using  the  grand  potential  and  NIM  for¬ 
mulation.  The  advantage  of  the  formulation  is 
that  it  is  applicable  to  the  ternary  liquid  system 
almost  without  change  except  for  the  subscript  i 
or  j  taking  values  1,  2,  and  3,  corresponding  to 
the  A,  B,  and  C  components. 


Fig.  2.  Liquidus  curves  of  III-V  semiconductors.  Solid  curves 
are  theory,  calculated  independently  by  Stringfellow  and 
Greene  (1]  and  Kikuchi  in  this  paper.  Circles  are  experi¬ 
ments. 

To  make  sure  of  the  relations,  and  also  to 
rewrite  equations  using  iv  and  fa,  we  now  re¬ 
examine  the  equations  in  section  3  one  by  one. 
Eq.  (3.1)  for  the  grand  potential  (j  holds  without 
change  when  is  for  i  *  1,  2,  and  3  and  as 
before  fa  is  a  parameter  related  to  the  chemical 
potential.  We  then  use  the  pair  approximation 
expression  of  the  CVM  and  write  u  as  in  (3.2). 
At  this  point,  however,  we  use  i„  introduced  in 
(3.13a)  and  fa  in  (3.15)  and  write  (5  as 

*  *  *  fro  X  ~  ^ 

2  *(*,)] 

x  [2  + 2  */M,] 

+0A[i-2y„].  (51) 
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Different  from  (3.13b),  we  need  three  energy 
parameters  *t2,  fa  and  <»• 

The  normalization  of  y*  is  written  as  in  (3.3), 
and  the  geometrical  relations  in  (3.4)  hold.  When 
the  grand  potential  is  minimized  with  respect  to 
y*  using  (3.4),  we  arrive  at 

In  (***)+«  Iny* 

-$0(A  +  A)-A0-O  (5.2a) 

or 

y*  -  eA*”CM/)°"-"'a”)  exp[-^  +  2“  (A  +  A)]  ■ 

(5.2b) 

The  set  of  equations  for  nine  y¥  values  is  solved 
with  ease  using  the  MIM.  The  NIM  is  particularly 
useful  when  the  number  of  variables  is  large. 

The  chemical  potential  n,  can  be  written  in  the 
same  way  as  in  (2.16).  Comparing  (5.2)  with 
(2.16)  we  see 

i- 1,2  and  3,  (5.3) 

which  are  similar  to  the  binary  case  (3.16). 

When  the  grand  potential  <3  in  (5.1)  is  a 
minimum  and  (5.2)  holds,  the  same  transfor¬ 
mation  as  (3.7a)  leads  to  the  identity 

A  -  <3/(2 N).  (5.4) 

Using  this  and  using  the  relation  between  <3  and 

F  in  (3.1),  we  can  derive 

F-< 3  +  (5.5) 

where  we  used  (5.3),  (5.4),  and  the  relation  be¬ 

tween  A  and  A  in  (3.15). 

Similar  to  the  last  comment  in  section  3,  we 
can  impose  one  condition  among  three  A  values 
because  of  the  normalization  relation 

Jjx,-1.  (5.6) 

It  is  convenient  to  impose  the  condition 

JjA-O,  (5.7) 

or 

Hs~  0.  (5.8) 


6.  Ternary  solid  phase 

As  an  example,  consider  a  ternary  solid  of 
composition  (GaAs)t-I(InAs)I.  This  possesses  a 
structure  modified  from  a  zinc-blende  structure, 
shown  in  fig.  3.  It  is  made  up  of  two  f.c.c. 
sublattice:  As  occupies  the  white  f.c.c.  sublattice 
in  fig.  3,  and  Ga  and  In  share  the  black  f.c.c. 
sublattice.  Let  us  write  the  three  species  Ga,  In, 
and  As  by  subscripts  i «  1,  2,  and  3  in  this  order. 

First,  we  examine  the  energy  of  this  system. 
Since  each  Ga  and  each  In  are  individually  sur¬ 
rounded  by  four  As  atoms,  the  nearest  neighbor 
interactions  contribute  the  energy 

E  -  S( 4«ijX,  +  4e2Jx2),  (6.1) 

where  N  is  the  number  of  black  sublattice  points 
in  fig.  3,  Xi  is  the  fraction  of  Ga  and  In  atoms  on 
this  sublattice,  and  el3  is  the  interaction  potential 
between  the  rth  and  the  3rd  (As)  species. 

Our  first  calculation,  based  on  the  energy 
expression  (6.1),  was  not  in  good  agreement  with 
experiment.  Therefore,  as  the  next  step,  we  have 
introduced  the  second-neighbor  interactions:  in¬ 
teractions  among  Ga-Ga,  Ga-In.  and  In-In  on 
the  black  sublattice  of  fig.  3.  We  write  the  frac¬ 
tion  of  i-j  (*',  /  *  1  and  2)  second-neighbor  pairs 
as  yv  and  the  interaction  potential  for  this  pair  as 
ty  Then  we  can  write  the  total  potential  energy 
as 

E  -  4Af  «,*  +  6N  £  jjj  <*y*  (6.2) 


Fig.  3.  Structure  of  lit— V  semiconductor  solid  phase.  If  all 
black  balls  were  the  same,  the  structure  is  the  rinc-blende 
type. 
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where  6N  is  the  number  of  second-neighbor 
bonds  in  the  black  sublattice  of  fig.  3.  Using  a 
procedure  similar  to  (2.4),  eq.  (6.2)  can  be 
rewritten  as 

E  =  N[4«,3  +  6<u)x1  +  (4*23  +  6*n)xi  +  24e,y12], 

(6.3) 

where,  analogous  to  (2.5),  we  define, 

4«,  =*  2«u  —  <n  —  <22-  (6.4) 

Since  we  assume  in  this  paper  that  no  vacancies 
occur  in  the  lattice  and  thus  that  the  white  lattice 
points  in  fig.  3  are  all  occupied  by  As  atoms,  we 
can  *«lte  die  statistics  of  the  black  sublattice 
only.  The  entropy  for  this  sublattice  can  be  writ¬ 
ten  exactly  as  in  (2.6)  with  the  coordination 
number  2at  ■  12.  When  F  ■  E  -  TS  is  minimized 
with  respect  to  yi2,  keeping  xt  and  x2  fixed,  we 
obtain  an  equation  with  the  form  of  (2.8): 

yuyi2=*yt2exp(4/3<,),  (6.5) 

the  only  difference  being  that  the  -e  in  (2.8)  is 
replaced  by  e,  in  (6.5).  The  definition  of  e,  in 
(6.4)  has  a  sign  different  from  (2.5).  Eq.  (6.5) 
can  be  solved  as  in  (2.11)  and  (2.12)  as 


r  ■  (1  +  4xIx2[exp(40<1)  -  1]}W.  (6.6b) 

When  the  free  energy  is  minimized  with  res¬ 
pect  to  yt2  and  (6.6)  holds,  the  chemical  poten¬ 
tials  of  the  species  i  =*  1  and  2  in  the  solid  phase 
are  written  in  exactly  the  same  way  as  in  (2.16). 
Comparing  the  energy  expressions  (2.4b)  and 
(6.3),  and  remembering  that  the  number  of  lat¬ 
tice  points  in  the  liquid  phase  in  section  2  was 
written  as  2 N  and  that  the  number  of  points  on 
the  black  sublattice  in  fig.  3  is  N,  we  see  that 

Mw  *  4*13  +  6e*  +  kT(6  In  y*  - 11  In  x(), 

(6.7) 

for  i  - 1  and  2, 

where  the  subscript  s  stands  for  the  solid  phase. 
At  this  point  note  that  the  free  energy  of  the 
black  sublattice  for  the  composition  xi  =  1  -  x2  is 
written  in  the  same  way  as  (3.10): 


In  the  special  case  where  X|-»l,yu  also  ap¬ 
proaches  unity  so  that  (6.7)  and  (6.8)  reduce  to 

FJ?/N-mS?*4<I3  +  6<i1,  when  X|-»l, 

(6.9a) 

and  similarly 

F«2/N  -  -  4<23  +  6<22,  When  x2  - 1.  (6.9b) 

When  we  examine  the  meaning  of  Ffl  we  see 
that  it  is  the  free  energy  of  the  GaAs  solid  phase 
at  the  temperature  of  interest  to  us.  Recall  that  it 
was  already  calculated  in  (4.4)  as  F£)^.  There¬ 
fore,  when  we  have  the  data  for  the  entropy  of 
melting  ASm  and  the  temperature  of  melting  Tm 
for  the  solid-phase  GaAs,  we  can  calculate 
F$/N  =  The  corresponding  data  for  InAs 
leads  to  the  value  Ffl/N  =*  Using  (4.4)  and 
(4.5),  we  can  write  them  explicitly  as 

mS-mIP+mI !!P-(TJ?>-T)432.  (6.10) 

Avf 

where  the  subscript  1  indicates  that  the  quantity 
is  in  the  liquid  phase,  and  /xf.V”  is  the  chemical 
potential  of  the  ith  species  (i  =  1  and  2)  in  the 
supercooled  liquid  state  in  which  As  composition 
is  50%  and  the  other  50%  is  the  ith  species. 

When  we  compare  our  (6.5)  with  eq.  (3a)  of 
S-G  [1],  we  see  that  our  4*,  in  (6.4)  corresponds 
to  S-G’s  fin c-ac  in  the  solid  phase  by 

4<,  =  2/Iac-sc/z,  (611) 

where  z  **  12  is  the  coordination  number  of  the 
f.c.c. 

7.  Geometry  of  the  ternary  coexistence  condition 

We  now  consider  the  condition  that  the  solid 
phase  of  section  6  and  the  liquid  phase  of  section 
5  can  coexist.  As  in  section  5,  we  write  Ga,  In. 
and  As  by  the  subscript  i  =  1,  2,  and  3.  When  the 
number  of  these  atoms  are  written  as  N>  (/  =  1,  2, 
and  3),  the  condition  in  the  solid  phase  is 

Ni  +  *2-N,-N,  (7.1) 

where  N  is  the  number  of  lattice  sites  in  one 
f.c.c.  sublattice.  The  total  number  of  atoms  in  the 
solid  phase  is  2N.  Thus,  for  the  liquid  phase  we 
require 


S"-; 


so 
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so  that  the  number  of  atoms  does  not  change 
when  the  solid  melts.  However,  relation  (7.1) 
does  not  hold  in  general  in  the  liquid  phase 
because  the  composition  in  the  liquid  phase  can 
be  different  from  that  in  the  solid  phase. 

Because  N  in  (7.2)  is  fixed  in  our  paper,  the 
number  of  independent  variables  for  the  com¬ 
ponents  is  equal  to  two  in  the  liquid  phase:  N\ 
and  N2.  The  allowed  range  of  Aft  and  N2  is  the 
rectangular  triangle  (0,0),  (2 N,0),  and  (0,2 N), 
shown  in  fig.  4.  For  the  purpose  of  our  analysis  it 
is  convenient  to  introduce  the  following  vari¬ 
ables: 

J/.N.  +  Ni;  V-iV.-Nj,  (7.3) 

or 

N.  -  *U  +  V);  N,  -  KU  -  V).  (7.4) 

In  the  solid  phase,  U  *  N  is  constant,  as  is 
shown  in  (7.1)  and  also  in  fig.  4,  and  thus  V  is  the 
only  composition  variable.  When  we  consider  the 
free-energy  axis  perpendicular  to  the  plane  of  N, 
and  N2  in  fig.  4,  the  free  energy  F,  for  the  solid 
phase  is  a  planar  curve  which  lies  within  the 
vertical  plane  U  ■  IV,  +  N2 »  N.  The  derivative 
of  F,  with  respect  to  V  is 

<7-5) 


Fig.  4.  Relation  between  the  W,  N{)  coordinate  system  and 
the  (l/,  V)  coordinate  system. 


where  we  have  used  the  for  the  solid  defined 
in  (6.7). 

At  the  composition  V"01,  the  free  energy  value 
is  written  as  F[0).  When  we  consider  a  plane 
which  contains  the  tangent  dFJdV  at  this  point, 
the  equation  of  the  tangent  plane  is  written  using 
an  arbitrary  constant  a  as 

F  -  F[0)  =  a(U  -N)+  (^)<0>(  V -  V<°>).  (7.6) 

This  arbitrary  coefficient  a  has  an  origin  similar 
to  that  of  the  a  used  in  (4.6)  for  the  binary 
case. 

Now  we  consider  the  free  energy  F|  surface  for 
the  liquid  phase  in  the  same  three-dimensional 
F-Ni~N2  space  based  on  fig.  4.  We  consider  a 
certain  composition  P  in  the  liquid  phase  as  in 
fig.  4  and  form  a  plane  tangent  to  the  F,  surface 
at  P: 

'-'r-wy-n 

+  ($l(V-V(P»-  C7-7) 

Note  that  in  this  equation  F  is  a  function  of  the 
variables  U  and  V  and  represents  a  plane  in 
F-U-V  space.  The  same  comment  holds  for 
(7.6). 

The  condition  that  the  solid  state  at  (N,  Vr,c”) 
and  the  liquid  state  at  (Uir\  V^)  coexist  is  that 
the  tangent  planes  in  (7.6)  and  (7.7)  become 
identical.  For  the  purpose,  two  sets  of  mathema¬ 
tical  equations  need  to  be  satisfied.  One  is  that 
the  coefficients  of  U  and  those  of  V  be  equal. 


respectively: 

/dFi\V) 

01  \dUJy' 

(7.8a) 

/d£Y0>  (9F |\« 

VdW  \*Wu 

(7.8b) 

The  other  requirement  is  that  (7.6)  pass  through 
the  point  (Fff>,  U(r),  Vn,))  or,  equivalently,  that 
(7.7)  pass  through  (F[°\  N,  V*0’). 

For  the  purpose  of  satisfying  these  two  con¬ 
ditions,  it  is  convenient  to  introduce  the  follow- 


R.  Kikuchi!  Theory  of  ternary  HI-  V  semiconductor  phase  diagrams 


51 


ing  two  functions: 

(7.9.) 

Then  we  can  prove  that  when  (7.8a)  and  (7.8b) 
are  satisfied, 

(7.10) 

satisfies  the  second  requirement.  The  proof  is 
similar  to  what  we  did  in  section  4.  Using  (7.9a), 
we  can  rewrite  (7.6)  as 


and  GP  of  (7.9)  we  first  examine  a.  We  use 
(7.14)  and  (7.5)  in  (7.8)  to  derive 

2Mn,  aisa) 

M>~mS  (7.15b) 

It  is  convenient  to  define  a  in  place  of  a  as 

«-«-M,+mS)).  (7.16) 

By  replacing  a  in  (7.15a)  by  a  and  adding  and 
subtracting  (7.15a)  and  (7.15b),  we  obtain 

(7.17a) 

a+M^-Mr-Mf-  (7.17b) 


<5‘°>  *  F  -  aU  -  (^)'°V.  (7.11a) 

When  we  use  (7.8),  we  can  rewrite  this  further  as 

,7ub) 

On  the  other  hand,  (7.10)  requires  that  (7.11b)  be 
equal  to  (7.9b).  Then,  when 

U  -  Uln  and  V  =  Vw,  (7.12a) 

it  follows  that 

F  =  FP.  (7.12b) 

This  means  that  (7.6)  passes  through  the  point 
(FT\  U°’),  V'r\  and  the  proof  is  finished. 

Thus,  what  we  need  when  calculating  the 
coexisting  solid  and  liquid  phases  is  to  evaluate 
<5(,0)  and  6P  in  (7.9).  The  derivative  dFJd  V  in 
(7.9a)  has  already  been  derived  in  (7.5).  The 
derivatives  in  (7.9b)  are  derived  using 

and 

(7.13) 


and  referring  to  (7.4)  as 


(7.14b) 

Before  substituting  these  derivatives  in  <5(,0) 


These  relations  are  helpful  when  we  examine  the 
limit  of  either  or  N2  approaching  zero.  For 
this  purpose,  we  first  rewrite  the  in  (6.7)  for 
the  solid  phase  as 

Mu  =  4«,3  +  6<a  +  fc7[ln  x,  +  6  In  (yjxf)], 

i  =  lor  2.  (7.18) 

When  the  solid  composition  x,  goes  to  zero,  we 
know  that 


and  therefore  that  when  Ni  approaches  zero.  Mu 
becomes  negative  and  its  absolute  value  becomes 
large  as 

HU-*kT  InXj.  (7.20) 

Since  we  expect  a  and  pP  in  (7.17)  to  remain 
finite,  we  see 

pP^kT lnxP  as  x<°>-0.  (7.21) 

The  relations  in  (7.20)  and  (7.21)  guarantees  that 
the  ternary  system  results  reduce  smoothly  to 
those  of  the  binary  system. 

We  now  rewrite  G[0)  in  (7.9a).  We  use  a  of 
(7.16)  and  dFJdV  of  (7.5)  in  (7.9a)  to  derive 

■  FJ°) —(<*+■ HmS> + mSSOIN 

-KmL<?)-mI3))V<o>  (7.22a) 

We  further  use  (7.1)  for  S  and  (7.3)  for  V*01  to 

transform 

(5i0)  =  F<°>  -6N-  (mS?W|0)  +  pWNW). 

(7.22b) 


37 


52 


R.  Kikuchil  Theory  of  ternary  IU-V  semiconductor  phcst  diagrams 


Since  we  see  from  (6.8)  that  the  terms  in  paren¬ 
theses  are  equal  to  F(Q),  we  arrive  at  the  simple 
result 

<5‘0)*-o N.  (7.22c) 

The  <5r’  in  (7.9b)  is  transformed  as  follows. 
We  use  (7.14)  for  the  derivatives  and  use  (7.3)  for 
U  and  V  to  derive 

6r  -  fv  -  oxr  -  m  wr  -  o*r  - 

(7.23a) 

Then,  using  the  expression  for  F\P)  in  (5.5),  we 
arrive  at  the  result 

GP-ltf'N.  (7.23b) 

Remember  that  IN  is  the  number  of  lattice 
points  in  the  lattice  model  of  the  liquid  phase. 

When  we  substitute  <5(0)  of  (7.22c)  and  of 
(7.23b)  in  (7.10),  we  obtain  the  relation  to  be 
satisfied  for  coexistence: 

-a=2tf\  (7.24) 

This  relation  further  leads  to  an  important  con¬ 
sequence.  When  we  use  this  in  (7. 17),  we  obtain 

M^-M^M^, 

+  (7.25) 

These  are  exactly  the  thermodynamic  conditions 
that  show  the  chemical  potential  of  the  AC 
molecule  (m2?0  in  the  solid  phase  is  equal  to  the 
sum  of  the  chemical  potentials  of  the  A  atom 
and  C  atom  in  the  liquid  phase.  The  fact  that  the 
coexistence  condition  (7.24)  leads  to  the 
established  thermodynamic  relations  (7.25)  sup¬ 
ports  the  correctness  of  the  analysis  in  these 
sections. 

8.  Derivation  of  ternary  liquid  us  and  solidus 

In  the  previous  section  we  formulated  using 
M’s-  As  discussed  in  section  3,  in  the  actual 
computation  we  use  pi's,  which  differ  from  fi’s 
by  A  +  otta,  as  was  shown  in  (3.16)  and  (5.3). 
Therefore,  the  next  step  is  to  rewrite  the  key 
equations  in  terms  of  p’s. 

Using  (5.3)  we  rewrite  (7.17)  as 

a  +  M(.?’  =  P?'  ~  PV  +  *>(«,v  -  <jj). 

for  i  =  1  and  2.  (8.1) 


We  can  let  pi  satisfy  either  (5.7)  or  (5.8).  In  this 
section  we  chose  the  latter.  Then,  (8.1)  can  be 
written  as 

p<lP)  -  a  +  nW  -  to(eu  -  e33),  for  i  =  1  and  2, 

PP  =  0.  (8.2) 

From  (6.7)  and  (6.9),  we  can  write 

Mi.i  =  m2?  +  kT(  6  In  y «  —  1 1  In  x, ), 

for  i  =  1  and  2.  (8-3) 

The  constant  term  m2?  is  obtained  by  using 
Vieland’s  formula  [5]  from  the  liquid  phase  of 
the  50-50  composition  as  in  (6.10).  The  quantity 
jtft3'  on  the  right-hand  side  of  (6.10)  is  the  chem¬ 
ical  potential  in  a  binary  liquid  phase,  so  that  it  is 
rewritten  using  (3.16)  as 

Ml.)3’  +  m(!53)  =  mI.'j)  +  A  +  tot,,  +  pk})  +  A  +  a>ei3. 

(8.4a) 

Since  we  use  (3.17)  to  calculate  p.2?,  we  can 
simplify  (8.4a)  as 

MtS3’  +  Mft”  =  2A  +  <o(c,,  +  e„).  (8.4b) 

This  is  a  part  of  m2?  in  (6.10).  Suppose  we  define 
Pi0/  as 

m2?  *  m2?  -  “(e„  +  e33),  i  =  1  and  2.  (8.5) 

Then  pi0/  can,  without  knowing  <=,,  +■  e33,  be 
evaluated  as 

A  C(i3) 

M2?  =  2A-(rS?)-7')=B-,  (8.6) 

rv  Avg 

where  A  is  calculated  for  the  liquid  phase  of  the 
50-50  composition.  This  relation  replaces  (6.10). 
Corresponding  to  (8.5),  we  define 

Pu  ■  M».i  -  «(e»  +  e33),  i  =  1  and  2.  (8.7) 

Then  the  relation  (8.3)  changes  into 

Pu  •  pf)+kT{6  In  yti  -  11  In  *,),  i  =  1  and  2. 

(8.8) 

We  then  use  (8.7)  and  rewrite  (8.2)  in  terms  of 

PW: 

P\f)  *  P(S)  +  a, 

M^)  =  m(,<?)  +  «,  (8.9) 

pr= 0, 
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where 


a  m  a  +  2w«33.  (8. 10) 

The  coexistence  condition  (7.24)  can  be  re¬ 
written,  vising  a  in  (8.10)  and  in  (3.16),  as 

—a  =2A.  (8.11) 

Note  that  20x33  cancels  out,  and  =  0  in  (8.9) 
is  used  here. 

The  long  transformations  done  so  far  in  this 
section  simply  mean  that,  for  the  phase  diagram 
calculation,  we  can  treat  the  liquid  phase  as 
though  the  A-A,  B-B,  and  C-C  pair  inter¬ 
actions  vanished. 

Thus,  the  actual  computational  steps  are  as 
follows. 

(a)  We  first  solve  1-3  and  2-3  binary  liquid 
equilibrium  states  for  50-50  composition  and 
then  estimate  (8.6): 

AS<i3> 

/2£,  =  2A-(n3>-7')^p-,  i  *  1  and  2, 

•<Va«i 

(8.12) 


where  A  is  the  Lagrange  multiplier  used  in  (3.6). 
As  was  the  case  in  section  4,  TS”  is  the  melting 
temperature  of  the  AC  solid,  and  dS^3'  is  the 
entropy  of  melting  of  one  mole  of  the  solid. 

(b)  Choose  a  composition  x\0)  -  1  -  x[0)  in  the 
solid  phase  and  calculate  (8.8): 

=  kT(6  In  y<P>  - 1 1  In  x<°>), 

i  =  1  and  2,  (8.13) 

where  yi°*  is  calculated  analytically  using  x\0)  as 
in  (6.6)  when  the  second-neighbor  interaction 
energy  e,  is  given. 

(c)  Assume  a  value  a  and  calculate  /Ip  from 
(8.9): 

m  ji<Q)  +  a,  for  i  =  1  and  2, 


(d)  When  ( i  =  1,  2  and  3)  are  thus  given, 

we  can  calculate  the  ternary  liquid  phase  using 
the  NIM  of  section  5.  The  normalization 
parameter  A  is  then  derived  and  used  to  see  that 
the  coexistence  relation  (8.11)  holds: 

—ft  =  2A.  (8.151 


Note  that  this  A  is  different  from  A  in  (8.12). 

If  the  relation  does  not  hold,  we  change  a  and 
repeat  the  procedure  from  (c).  When  (8.15) 
holds,  then  the  liquid  phase  for  this  a  is  the 
phase  which  coexists  with  the  solid  phase  chosen 
in  (b). 

As  an  example,  we  calculate  the  phase 
diagram  for  Ga-In-As.  The  calculation  is  based 
on  the  binary  phase  diagrams  in  fig.  2.  which  were 
calculated  using  the  parameters  in  table  I.  The 
additional  parameter  we  need  for  the  ternary 
case  is  <,  in  (6.4)  for  the  second-neighbor  inter¬ 
action  in  the  solid  phase.  In  view  of  (6.11),  we 
used  the  value 

4e,  =  2800/6  cal/mole  (8.16) 

which  is  consistent  with  the  value  used  by  S-G. 
/2(GaAs-InAs)  =  2800  cal/mole.  Our  es  cor¬ 
responds  to  S-G’s  fi  divided  by  24.  The  positive 
value  of  e,  in  (8.16)  shows  from  (6.4)  that  Ga  and 
In  in  the  solid  have  a  tendency  to  repel  each 
other.  A  liquidus  isotherm  at  1250  K  and  a  few 
tie  lines  are  plotted  in  fig.  5.  The  solid  phase  is 
always  on  the  pseudobinary  line 
(GaAs)i-,  (InAs),. 

The  liquidus  and  solidus  on  the  pseudobinary 
line  are  plotted  in  fig.  6a.  The  calculated  curves 
are  exactly  the  same  as  those  calculated  by  S-G 


Fig.  5.  The  isothermal  liquidus  of  Ga-In-As  ternary  system 
at  1250  K,  calculated  independently  by  Stringfellow-Greene 
[1]  and  Kikuchi  in  this  paper.  The  numbers  along  the  curve 
are  x  in  the  solid  composition  (GaAs)i-,(InAs),  Four  tie 
lines  arc  shown 
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Fig.  6.  Pseudobinary  liquid  vs  and  solidus  curves  of  InAs- 
GaAa  (a)  and  InSb-InAs  (b).  The  solid  curves  were  cal¬ 
culated  independently  by  Stringfellow-Greene  [1]  and 
Kikuchi  in  tins  paper,  and  the  cr.des  and  triangles  represent 
experimental  results. 


and  agree  well  with  the  experimental  points 
marked  in  the  figure. 

Although  the  computed  results  are  the  same,  a 
slight  difference  in  the  identification  of  the 
energy  parameter  «,  between  this  paper  and  S-G 
may  be  worth  pointing  out.  We  identified  «,  as 
the  parameter  for  the  second-neighbor  inter¬ 
action  (in  the  zinc-blende-type  lattice)  between 
Ga  and  In  in  the  solid  phase,  as  introduce*,  in 
section  6  and  illustrated  in  fig.  3.  On  the  oth'r 
hand,  S-G  regarded  /l(GaAs-InAs)  as  the 
parameter  for  the  GaAs-InAs  interaction. 

The  calculation  was  also  done  for  the  In- As 
Sb  system.  The  energy  parameters  are  those  in 
table  I  together  with  the  solid  phase  value 

4«,-  2900/6  cal/mole.  (8.17) 

The  liquidus  and  solidus  on  the  pseudobinary 
line  are  plotted  in  fig.  6b;  they  are  in  good 
agreement  with  S-G’s  curves  and  with  experi¬ 
mental  data. 


9,  The  orthogonality  relation 

At  the  end  of  section  S  we  commented  that 
either  (5.7)  or  (5.8)  could  be  used  as  the  con¬ 
dition  for  the  fi..  In  section  5,  we  used  /I3  =  0  in 
(5.8).  We  now  form  linear  combination  of  fi,  and 
define 


0*i  +  Mi)/3,  i=  1,2  and  3,  (9.1) 

so  that 


which  is  of  the  form  of  (5.7).  Solving  (9.1)  for  /2,'s 
we  can  write 


fin  =  /*?  - 45,  f- 1,2  and  3.  (9.3) 

Since  the  /Z*  satisfy  condition  (9.2),  we  can 
plot  them  on  a  star-shaped  graph.  An  example  is 
shown  in  fig.  7.  Each  /Z*  (i  =  1,  2  or  3)  point 
corresponds  to  a  coexisting  liquid  state  P  and  a 
solid  state  Q  and  thus  to  the  tie  line  PQ  We  can 


Fig  7  The  chemical  potential  diagram  corresponding  to  fig  5 
The  numbers  along  he  curve  are  the  »  values  corresponding 
to  those  n  ftg .  5  The  m  curve  it  x  is  perpendicular  to  the  tie 
line  at  x  in  fig  5 
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now  prove  that  the  tie  line  PQ  is  perpendicular 
to  the  tangent  of  the  /Z*  curve. 

From  (5.3)  the  chemical  potential  p.f)  at  the 
liquid  point  P  can  be  written  as 

=  jZj  +  A+<u€jj  =  /i*-/Z$  +  A+  wea, 

i  *  1  2  and  3.  (9.4) 

where  we  used  (9.3).  The  chemical  potential  for 
an  AC  unit  and  a  BC  unit  in  the  solid  phase  can 
be  written  as  pff  and  /xSS*  in  (7.25).  Furthermore, 
using  (9.4)  we  can  write  them  as 

pi?  =  /!! -/2*  +-2A  +  <u(«a  +  <33),  i  =  1  and  2. 

(9.5) 

Along  a  liquidus  isotherm  for  T,  the  change  of 
the  grand  potential  <3  due  to  the  change  in  the 
chemical  potential  is  written  for  the  liquid  phase 
as 

Aar=2N^xr^r,  (9.6a) 

where  21V  is  the  total  number  of  atoms  in  the 
liquid  phase.  Using  (9.4),  we  can  rewrite  this  as 

d<3r>  =  2 N  [Jj  xT'AH:  -  4AS  +  4A  ].  (9.6b) 

At  the  other  end  of  the  tie  line  the  solid  phase 
is  made  up  of  Nxu  AC  units  and  NxtJ  BC  units; 
these  have  chemical  potentials  jtjfP  and  tx($\ 
respectively,  as  shown  in  (9.5).  Thus,  the  change 
of  the  grand  potential  along  the  solid  phase 
isotherm  is 

=  (9.7a) 

In  rewriting  this  expression,  note  that  the  coor¬ 
dinate  x(°>  for  the  solid  state  in  the  Gibbs  trian¬ 
gle  plot  is: 

*i0)  *  i  -  1  and  2, 

*P-l  (9.7b) 

Combining  these  relations  with  (9.5),  we  can 

rewrite  (9.7a)  as 

d(5<°>  -  2N[x\<»Aji  T  -  A/l  5  +  A\  ].  (9.7c) 

The  grand  potentials  <j\P)  and  <3<0>  are  equal 
when  P  and  Q  are  connected  by  a  tie  line.  Thus, 


A<j\P)  and  A<j{P)  are  also  equal.  Then,  by  equat¬ 
ing  (9.6b)  and  (9.7c),  we  obtain 

^(xf>-x<°>)dM!=0.  (9.8) 

Note  that  xf*  —  xi0)  and  Afi.*  are  components  on 
three  coplanar  axes.  It  was  proved  in  the  appendix 
of  ref.  4  by  simple  geometry  that  the  components 
of  these  vectors  on  Cartesian  coordinate  systems 
obey  an  equation  similar  to  (9.8).  Therefore,  we 
can  conclude  that  the  tie  line  represented  by 
x^-xi01  (»  =  1,  2,  and  3)  is  perpendicular  to  the 
tangent  of  the  (L  ?  curve. 

Appendix 

Proof  that  the  grand  potential  decreases  as  the 
interaction  proceeds  in  the  SIM 
In  this  appendix  we  prove  that  the  grand 
potential  <3  or  <P  in  (3.2)  decreases  in  one  itera¬ 
tion  step  of  the  NIM.  We  rewrite  (3.5)  as 
i Ity  ■  /3wev  —  j(2a»  —  1)  In  (xpc,)  +  at  In  y, 

-&(*+&)-  A0=O,  (A.l) 

where  x,  and  x,  are  the  input  to  the  iteration 
cycle,  %  is  the  output,  and  other  quantities  are 
constants.  We  can  rewrite  the  grand  potential  <t> 
in  (3.2)  as 

0{y«}  =  2  evy  *-&*>- 1) 

**/ 

x  ^2  x<  In  x,  +  2  x,  In  x,] 

+«.  2  y#  *"  y « -  \P  [2  *& +Xv,]. 

(A. 2) 

This  is  a  function  of  the  input  x,  and  y,,  values,  as 
indicated  by  {y„}.  Note  that  in  (A.2)  x,  and  x, 
represent 

Xi-^yu  and  x,  =  £  y„,  (A.3) 

In  later  transformations  we  use  <P{y„},  which  is 
derived  by  replacing  the  x ,  and  y„  values  by  x, 
and  %  values.  The  relations  (A.3)  hold  for 
quantities  with  carets  also. 

Next,  we  use  (A.l)  and  form 
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x  [2  x< ,n  x< + 2  x< ln  x<  ] 

+«  2  y<i |n h -te  [2  x&  +  2 *#/] 

-A  0-0.  (A.4) 

Since  this  vanishes,  we  can  simplify  <t>{ya}  in 

(A.2)  by  subtracting  (A.4)  from  (A.2)  as 

<0{y,}  *  «  2  y»0n  y*  -  in  9a) + a/s.  (a.5) 

V 

Then  from  (A.l)  we  form 

2  9a+v  2  *»9a  ~  i(2<o  - 1) 

W  Q 

x  A  In  Xj  +  2  2/  In  x,] 

2  9* ,n  9a~  [2  +  2 

-A0 -0  (A.6) 

and  subtract  this  from  <P{yv},  which  is  derived 
from  (A.2): 

*{*,}« +2w-l) 

x  ^  i,(ln  i,  -  In  x.)  +  2  2/ On  2,  -  In  x,)j  +  A0. 

(A.  7) 

We  note  that  ♦{y#}  in  (A.5)  is  the  input  value 
of  the  grand  potential,  and  ♦{$,}  in  (A.7)  is  the 
output  value.  We  are  interested  in  the  difference 
between  the  two: 


♦{y*}  -  {/*}  =  (2<o  -  1)  2  *.  ln(x, — x() 

I 

+«  2  ln(y#/y*)-  (a.8) 

In  transforming  this  expression,  we  use  Gibbs’s 
lemma  [7]: 

dTje^'  +  l-e^aO,  (A.9) 

which  holds  for  any  real  number  Art.  The  left- 
hand  side  is  a  decreasing  function  of  Art  for 
negative  values  of  Art  and  an  increasing  function 
of  Art  for  positive  values  of  Art  ;  it  vanishes  for 
Art  =  0.  We  rewrite  (A.8)  as 

*{y#}  “  *  (2w  - 1)2  ^  ln(Xj/x()  +  x(  -  xf] 

I 

+  "  2  [y«  toOVy.)) + %  -  ya\ ■  (A- 10) 

V 

Since  the  summand  of  each  term  is  non-negative 
because  of  (A.9),  we  therefore  have  proved  that 

4>{y#}-*{y#}&0.  (A.11) 
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ABSTRACT.  The  pair  approxlaation  (which  Is  equivalent  to  the  quasl-chemlcal  method)  Is  used 
to  calculate  llquldus  of  II-VI  compound  semiconductors.  For  the  liquid  phase,  a 
quasl-lattlca  structure  Is  assumed,  and  the  pressure  and  the  association  property 
are  taken  Into  account.  Vleland's  formula  Is  used  for  the  free  energy  difference 
between  liquid  and  solid.  Numerical  calculation  was  made  with  the  Hg-Te  and 
Cd-Te  systems.  For  each  system,  one  set  of  parameters  Is  found  which  computes 
the  llquldus  In  good  agreement  with  experiments  for  the  entire  composition  and 
temperature  ranges. 

A  particular  attention  la  paid  to  the  shape  of  the  llquldus  curve  near  its 
top.  Experlmsntally  It  looks  pointed,  although  thermo dynamic!  requires  It  to  be 
rounded.  Our  calculation  shows  that  when  the  AB-A  and  AB-B  Interactions  are 
repulsive,  the  theoretical  curve  looks  pointed  and  the  top  region  has  a  very 
smell  radius  of  curvature. 

1.  Introduction 

A  theory  of  llquldus  of  III-V  semiconductors  was  formulated  using  the  quasl-chemlcal 
method  by  Strlngfellow  and  Greene  (1),  and  later  by  the  present  author  (2).  In  the  latter, 
the  cluster  variation  method  (CVM)  approach  (3,4)  was  used;  the  pelr  approxlaation  of  CVM 
Is  equivalent  to  the  method  commonly  known  as  the  qussl-chnmlcal  method  (5) . 

In  the  present  paper,  we  extend  the  previous  trestment  (2)  to  the  II-VI  semiconductor 

llquldus  calculations.  The  mein  dlfferencs  Is  the  aseocletlon  property,  that  means 

Inclusion  of  molecular  species  In  the  liquid  phase  In  the  II-VI  case. 

The  liquid  phase  Is  represented  by  a  lattice  structure.  In  the  pair  approximation  we 

are  using  In  the  paper,  the  Information  on  the  structure  of  the  lattice  appears  only  In  the 
coordination  number  which  Is  written  es  2w. 


2.  Chemical  Potentials  of  the  Liquid  Phase 

The  two  atomic  species  In  the  system  are  written  as  A  and  B.  They  are,  for  example, 
Hg  and  Te.  He  assume  that  In  the  liquid  phase  the  molecular  species  AB  also  exists.  We 
use  a  lattice  model  for  the  liquid  phase,  and  these  species  can  sit  only  on  lattice  points. 
Although  In  reality  obviously  the  liquid  does  not  have  a  lattice  structure,  the  approxi¬ 
mation  Is  considered  eccepteble  because  the  local  correlation  among  atoms  and  molecules 
In  the  reel  liquid  system  can  be  simulated  adequately  using  the  lattice  structure,  as  far 
as  no  superlattice  structure  Is  Introduced. 

The  three  species  A,  B,  and  AB  are  designated  by  i-l ,  2,  and  3.  We  also  include 
vacancies  (t“4)  la  the  lattice  In  order  to  take  Into  account  the  pressure  effect.  The 
probability  of  finding  an  1th  species  on  a  lattice  point  Is  written  as  x.  and  the 
probability  of  finding  an  l-J  nearest-neighbor  pair  Is  written  as  y,,  These 
variables  satisfy  the  geometrical  relation  or  the  reduction  relation: 
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l-l,...,4 


Since  there  le  no  superlattice  structure  In  the  lattice,  the  symmetry  relation 
holds: 


*u  ■  *n 


The  normalization  relation  t« 


Before  going  Into  the  eeln  treatment  of  the  liquid  phase,  In  this  section  we  discuss 
the  chealcal  potentials  of  species  In  the  liquid  phase.  For  this  purpose  we  write  the 
Rmlmholts  free  energy  F»E-TS.  The  energy  E  Is  written  as 


'Y.  Y.  cl^tj  -  Me3*3 
l-l  J-l 


where  B  Is  the  number  of  lattice  points  In  the  liquid  systea,  ct4  Is  the 
Interaction  potential  for  the  nearest-neighbor  1-J  pair,  and  e3  Is  the  formation 
energy  of  the  AB  molecular  species.  Since  we  assume  no  interaction  energy  with 
vacancies,  we  require 


el*  *  C41  -  0  for  1-1 • 


The  entropy  foraula  Is  the  saaa  as  Equation  2.6  In  Reference  2: 


S  -  U»  (2wl)  E  l(*t)  -  »  l(ytJ) 


where  the  L(x)  function  originates  In  the  logarltla  of  a  factorial  and  Is  defined  as 

l(x)  *  x  In  x  -  x  (2.7) 

The  3  expression  In  Equation  2.6  la  for  the  pair  approximation  of  the  CVM  and  la  the  same 
as  the  expression  used  In  the  quasi -chemical  method  (3). 

When  we  use  E  In  Equation  2.4  and  S  In  Equation  2.6,  we  can  write  the  Helmholtz  free 

energy  F  as 


f  -  a  £  *11*11  "  *3*3  "  kT  *  “  1C  UTij) 


For  the  discussion  of  chealcal  potentials,  it  Is  convenient  to  use  the  relation  of 
Equation  2.1  and  write 


(2.-  -/ ; 
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yu  ■  *1  -  £  yu 

Bated  on  Chit  rtlacion,  we  Crete  f  at  a  function  of  xj  (1*1 ,...,*)  and  y4j 

(l*J>. 

The  run  bar  of  Che  1th  epaclaa  In  a  lattice  la  written  at  N1(  which  tatlsflea 


Nt  -  xt  N  ,  1-1,...,* 


(2.1 


Note  N,,  la  Che  nuaber  of  vacanclaa.  Whan  we  uaa  the  Helmholtz  free  energy  F  at  the 
theraodynaalc  potential,  the  total  voluaa  of  the  ayaeea  la  fixed.  In  the  lattice  nodal, 
thla  meant  that  the  total  nuaber  of  lattice  points  N  la  fixed.  Me  can  choose  Nt,  N2,  ai 
Nj  as  Independent  and  write  N4  at 


N4  -  S  -  Nt  -  N2  -  Nj 


(2.1 


Me  then  calculate  chemical  potentials  aa  follows: 


“I  5  Wn2.n3  ’  Wn2,n3,n4  •  Wn1,n2,n3 


(2.: 


When  we  use  the  free  energy  expreaalon  In  Equation  2.8,  we  can  obtain  explicitly, 


where 


ul  “  "*11  +  kT|ln  +  u  ln^rn/xi2^j  "  \ 

“v  3  I54  *  kT[ln  x4  +  “  ln(?**/x*2)j  • 


(2.: 


(2. 


We  can  perform  the  saae  differentiation  aa  Equation  2.12  for  1-2  and  3  and  obtain 


M1  "  “*11  ~  ci  +  bljln  xt  +  «  ln(7ti/*iZ)J  "  >V  for  i“l>  2  and  3  (2.: 

where  we  wrote  tf  aa  a  genaral  tarn  but  It  la  defined  aa 

*i  -  0  for  1  a  3  (2.16) 

and  c3  la  the  energy  paraaater  In  Equation  2.*  and  2.8b. 

The  total  nuaber  of  A  atoms  Including  those  In  molecules  AB  la  written  as  NA.  We 
define  the  corresponding  quantity  Ng  for  Che  B  atoaa.  Then 

NA  -  Nj  +  N3  , 

NB  -  N2  +  N3  .  (2.17) 

When  the  compositions  NA  and  Ng  are  given,  tha  system  adjusts  itself  Co  coma  to  the 
value  of  Nj  such  that  the  free  energy  la  a  minima  with  respect  to  N3.  In  order  to  see 
thla  relation,  we  can  use  N4>  SB>  and  N3  as  Independent  and  write  Nj,  N2,  and  N4  as 

N1  “  nA  *  N3  • 

N2  -  NB  -  i»3  , 

N4  -  N  -  N.  -  N.  +  N,  .  (2.1! 
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When  we  minimise  f  with  respect  to  Nj  heaping  SA  and  N*  fixed,  we  can  use 
Equation  2.12  and  corresponding  ralatlona  for  1-2  and  3  to  obtain 

U}  -  lit  +  u2  .  (2.19) 

This  la  tha  equilibrium  condition  aaong  atoaa  A  and  B  and  tha  molecular  spociaa  AB. 

Ha  can  uaa  Equation  2.18  again  and  derive  that  ur  In  Equation  2.12  la  identical  to 
tha  following  differentiation: 


“1 


(2.20) 


Thla  aaya  that  the  chemical . potential  u1  can  be  calculated  by  either  Bquatlona  2.12  or 
Equation  2.20. 


3 .  Grand  Potential 

In  nuaerlcally  solving  tha  liquid  state.  It  la  convenient  not  to  fix  tha  composition 
Hi  but  to  let  It  vary  by  using  Lag rang a  multipliers.  We  define  the  grand  potential  G 

aa 


2  -  E 


(3.1) 


where,  at  this  stage,  u|  is  regarded  aa  Lagrange  multipliers.  Since  N  Is  fixed  and  we 
treat  !t4  aa  dependant  as  In  Equation  2.11,  we  actually  do  not  need  u4  In  Equation  3.1. 
So  we  require 


“4 


0 


We  use  Equation  2.9  and  rewrite  the  energy  E  in  Equation  2.4  as 


irftere  we  define 


Special  eases  are 


and 


ei4 


'4i  ■  *  2  cii 


(3.2) 


(3.3) 


(3.4a) 


(3.4b) 


(3.4c) 


When  we  combine  E  in  Equation  3.3  and  Che  chemical  potential  terras  in  Equation  3.1, 
we  can  write 


—  14 

S  M-Ni  -  c\.  ytj  ♦  »  S  C«*ii  -  ‘i  -  (3.5) 


i-1 


t.J 


i-1 
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We  can  simplify  this  when  we  write 

UJ.  s  -  "ii  ♦  ei  • 


(3.6) 


Because  of  Equation*  3.2,  2.5,  and  2.16,  we  see 

S4  *  o 


(3.7) 


Uhan  we  combine  Equation*  3.5  and  2.6,  we  can  write  the  grand  potential  G  in 
Equation  3.1  explicitly  a* 


♦  5  £  -  Sw2Ze;iy; :  -  6  2  UjXj^ 

"  ij  lJ  l  1  l 

-  (2<a  -  1)  S  L(x,)  +  »  Z  L(y,:)  *  8Ml  "  S  y  L  j ) 
i  1,1  J  \  l.J  -7 


(3.8) 


where 


and  X  ia  a  Lagrange  multiplier.  The  X  term*  are  added  in  Equation  3.8  in  order  to  Cake 
into  account  the  normalixation  of  yjj . 

^  The  equilibria*  atate  i*  obtained  by  minimizing  G  in  Equation  3.8  with  T  and 
fixed: 


rij 


5  8*i  t"  -  ly  'l  In  (XjKj)  ♦  wlnyj^j  -  6  (»£  ♦  Sj)  -  X8  ”  0  (3.10) 


At  thia  point  we  digreaa  and  examine  the  i«j  caae.  Using  Equation  3.4  we  can  derive 
Sj  ♦  X  »  kT^ln  x^  *  uln (yii^xi2)}  for  *•  ”  i >2 ,3  .  (3.11) 

The  fourth  equation  can  be  simplified  by  uaing  Equation*  3.7  and  3.4  a* 

X  •  kr|lnx4  ♦  m  In  (y^/*^  2)j  .  (3.12) 


Hhen  we  compare  thia  with  Equation  2.14,  we  identify  u  uaed  in  Section  2  aa 


Uy  -  X  '  (3.13) 

Ua*  of  thia  in  Equation  2.15  and  compariaon  with  Equation  3.11  allow  us  to  derive 

-  “*u  ♦  ,  i  •  1,2,3  .  (3.14) 


Vdrthar,  when  we  co^ara  thia  with  Equation  3.6,  we  can  identify 

»i  ■  ,  i  -  1,2,3  .  (3.15) 

which  mean*  that  uaed  in  the  expreaaion  of  G  in  Equation  3.1  i*  actually  the  chemical 
potential  itaalf.  Therefore,  from  now  on  wa  can  drop  the  prime  from  u[. 

„Whan  thia  relation  ia  uaed,  £wt»i  become*  the  Gibb*  free  energy  G  itaelf. 

Thu*  G  in  Equation  3.1  i*  1 
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where  v  ia  Che  voluaa  par  lattice  point  ia  tha  lattice  nodal  of  che  liquid  phaee. 

Tha  identity  ia  Equation  1. 15  leads  to  another  inportent  relation.  Ue  can  write  the 
aeaociation  equilibria  relation  in  Equation  2.19  ia  tenea  of  u:  using  Equations  3.15  and 
3.14: 


»1  ♦  5j  -  S3  -  u  (e33  -  eu  -  t22)  •  (3  5u  (3.17) 

We  can  interpret  the  aeaning  of  tha  energy  constant  u  as  follows.  When  in 
Equation  3.11  is  written  as 


Si  s  IcT  In  at  ,  i  -  1,2,3  ,  (3.18) 

then  a^  ia  called  the  activity.  Equation  3.11  together  wit.h  3.12  gives  the  explicit  forte 

of  ai  as 


which  says  that  a{  ia  roughly  proportional  to  the  fractional  aaount  xf.  When  ue  use 
*i,  we  can  write  Equation  3.17  as 


■  asp  (-u/kT)  (3.20) 

This  shows  that  the  energy  constant  u  controls  che  aaount  of  association,  a  larger  -u 
leading  to  a  larger  a],  i.e.,  larger  association. 

When  •  in  Equation  3.8  is  a  ainiaua  and  Equation  3.10  holds,  ue  can  simplify  *  by 
foraing 


(3.21) 


this  weans  that  all  teras  in  *  cancal  each  other  except  the  only  one  81.  Coa^arison  of 
Equation  3.21  and  3.18  leads  to  che  relation 


1  "  -V  •  (3.22) 

How  we  go  back  to  the  equations  in  3.10,  and  write  it  as 


*ij  "  -  Siij  ♦  fj;  (1^  ♦  V]  (xix j)(2uHl)/(2a) 


Whan  8  ■  1/kT  and  Si  are  given,  chit  relation  and  the  reduction  relations  in  Equation  2.1 
fora  a  set  of  aiaultaneoua  equations  to  be  solved  for  yjj.  We  can  solve  thaa  using  the 
MZM(3) . 

The  energy  paraaeters  we  need  are  seven:  el2,  c,,,  e23,  el4,  t24,  t31,  in 
Equations  3.4  and  u  in  3.17.  The  values  of  these  seven  paraaeters  are  to  be  supplied  by 
whetever  naans  before  the  coaputation  starts.  In  solving  che  sec  of  siaulcsneous  equations 
we  fix  two  values  ut  and  fi2.  When  u,  and  S2  are  given,  u3  is  deterained  froa  the 
association  relation  Equation  3.17;  note  v„  ■  0  at  in  Equation  3.7.  We  then  proceed  to 
solve  the  equations  using  HIM.  After  tha  iteration  has  converged,  che  pressure  is  derived 
froa  Che  relation  in  Equation  3.22.  Tha  two  adjustable  paraaatert  Q.  and  u2  can  be  so 
chosen  that  the  pressure  p  and  che  coapoaicion  HA/Hg  can  take  chosen  values. 


4,  Derivation  of  the  liouidus 

Derivation  of  the  liquidus  curve  was  fornulated  in  Section  4  of  Reference  2  bated  on 
the  geoaetrieal  construction  thet  che  tangent  of  the  Gibbs  free  energy  curve  for  che  liquid 
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phase  goaa  through  cha  Gibbs  free  energy  point  for  the  solid  phase.  We  can  use  the  siitilar 
formulation  in  tha  praaent  problem.  Ue  assume  as  in  Raferenca  2,  that  the  solid  phase  has 
tha  fixed  composition,  Ng  *  Ng,  and  thus  the  Gibbs  free  energy  for  the  solid  at  T  is 
raprasaatad  by  a  single  poiat  rather  than  a  curve. 

Tha  geometrical  construction  of  Raferenca  2  leads  to  tha  coexistence  condition  that 


V"  ♦  «2U)  *  »1(,> 


(s) 


(4.1) 


share  the  superscripts  ( l)  and  (s)  stand  for  liquid  and  solid,  respectively.  This  equation 
holds  in  the  present  case  also. 

The  chemical  potentials  for  the  solid  state  can  be  calculated  using  Vieland's 
method  (6),  and  is  uritten  as 

„/•>  ♦  -  Uj(c)  ♦  Uj<c)  -  (T„  -  T)  t£m  (4.2) 


where  (c)  indicates  the  supercooled  liquid  state  of  the  composition  Ng  ■  Ng  ac  T,  T„ 
is  the  melting  temperature  of  the  solid  end  A.%,  is  the  entropy  of  melting  per  molecule 
Ah.  The  value  of  Tm  is  to  be  supplied  from  experiments. 

In  the  present  case  we  can  use  Equation  2.19  and  aiisply  write  u3  for  u t  *  uz  in 
Equation  4.2.  Further,  wa  can  use  Equation  3.14  to  write 


PjU)  -  -  ij(c) 


(4.3) 


Note  a  caret  on  u  on  the  right-hand  side.  Combining  these  chree  equations,  we  can  write 
the  condition  for  the  liquidua  aa 

ij<1)  -  H3<c)  -  (Ta  -  T)  as,,  .  (4.4) 


Thus  tha  procedure  of  deriving  the  liquidua  curve  ia  tha  following.  At  Che 
temperature  T  and  pressure  p,  the  a  iparcooted  state  (c)  ia  first  solved.  This  solution 
gives  the  value  S3'c^.  Then  Equation  4.4  is  used  to  derive  for  the 

liquid  state  which  coexists  with  the  solid  st  this  T  and  p.  Equation  3.17  is  then  used  to 

evaluate  •  Hh,n  «»*ign  *  value  to  -  u2'*J, 

all  the  parameters  ace  ready  to  solve  the  simultaneous  Equations  3.23  for  the  liquid  state 
(*).  The  adjuatabla  value  is  determined  so  chac  the  pressure  of 

the  state  becomes  equal  to  the  desired  value  p. 

5._  Choice  of  the  Energy  Parameters 

In  Section  3  we  noted  that  the  prasant  formulation  has  seven  energy  parameters  ac  our 

dispoaal.  Among  the  seven,  tha  three  tlt>  c2l|  and  c34  are  related  to  vacancies  and  hence 

have  little  effect  in  tha  liquidua  caleulationa  because  Che  number  of  vacancies  (in  the 
pseudo-lattice  model)  ia  always  small,  lass  chan  10"*,  in  the  present  calculations. 

In  the  III— V  case,  in  which  no  molecular  species  exiscs,  tj3  and  e23  do  not  appear  in 

tha  theory  and  hence  ct2  ia  the  only  energy  parameter.  Ha  know  in  the  III-V  case  (2)  that 

Cj2  controls  tha  shapa  of  the  liquidua  curve,  making  it  higher  or  lower.  In  the  present 
II-VI  case,  also,  cl2  controls  the  liquidua  curve  ia  the  similar  way.  Whan  -c12  is  large 

tha  liquidua  goes  up,  excapt  at  Tg  which  is  fixed. 

The  energies  c13  and  c2>,  together  with  u,  control  the  amount  of  association.  In 
order  to  single  out  the  association  property,  wa  calculate  the  special  case  c13a  c23, 
keeping  cha  liquidus  symmetric.  The  calculated  results  are  in  n«- i-  in  which  the  number 
written  next  Co  a  curve  ia  Che  value  of^che  ratio  i13/ cl2  ■  «2j/*i2.  For  lower  curves  for 
which  the  ratio  is  positive,  and  hence  i13  *  i23  is  negativs  (i.e.,  attractiva),  the 
liquidus  shows  the  rounded  top.  For  upper  curves  for  which  the  ratio  is  negstive,  i.e., 
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e13  «  e23  is  positive  and  repulsive,  Che  liquidus  shows  s  pointed  peek.  As  the  shepe 
changes  iron  the  rounded  top  to  the  pointed  top,  around  -2  of  the  ratio,  the  shoulders  are 
si no at  linear.  At  -3  and  -4  of  the  ratio,  the  inflection  points  are  observed.  For  -5 ,  a 
phaee  separation  occurs  in  the  liquid  phase.  The  liquidus  similar  as  the  -5  case  is 
observed  in  the  Zn-rich  side  of  the  2n-Te  binary  liquid  (7). 

Although  the  upper  curves  in  Fig.  1  are  pointed,  the  peak  still  has  a  rounded  top 
with  a  small  radius  of  curvature.  This  is  in  agreement  with  the  thermodynamic 
requirement  (8) . 


Hg-Te  binary  liquidus.  Dots 
are  experiments  (9). 


Liquidus  curves  for  different 
values  of  atom-molecule  Inter¬ 
action  parameters.  A  horizon¬ 
tal  short  line  narks  600*C  for 
each  curve. 


Replots  of  the  low  temperature 
parts  of  the  Hg-Te  liquidus. 
The  vertical  axis  le  the 
composition  of  the  minority 

species. 
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Figure  2  shove  Che  Hg-Te  binary  liquidut.  The  does  are  Harman's  experiments  (9).  In 
calculating  the  theory  vs  used  Che  following  values.  The  melting  temperature  is  Harman's 
experimental  value  1,  *  673*C.  The  entropy  of  melting  at  Tn  is  ASq  ■  9.2  e.u.  which 
is  in  Laugier's  paper  (10).  ^Other  energy  values  are  £12“  -200/6  cal,  £13/612  "  -8, 

*23/632  *  —0.5,  *i4  m  *24  *  *3%  m  "1000**  and  u  ■  —3000  • 

Figure  3  shows  the  low  temperature  parts  of  the  liquidus  curves  for  the  Hg-rich  side 
and  the  Te  side.  The  Te  values  on  the  Hg  side  are  based' on  Riley's  experiment  (11).  They 
are  plotcad  in  the  logarithmic  scale  so  that  a  small  asiounc  of  the  minority  component  is 
seen  more  accurately  than  in  Fig.  2.  The  curves  in  Fig.  3  are  close  to  linear.  It  is 
found  by  crying  many  parameter  value  combinations  that  the  slope  of  the  line  in  Fig.  3  is 
controlled  by  ASn  and  almost  solely  by  it.  For  a  certain  value  of  ASm,  when  we  change 
other  energy  parameters,  the  line  in  Fig.  3  goes  up  or  down  but  the  slope  remains  almost 
unchanged  at  Che  value  determined  by  ASq.  The  docs  in  Fig.  3  are  experiments  and  the 
curves  are  the  theory;  they  agree  well  and  the  theory  passes  Che  test  which  is  more  severe 
chan  Fig.  2.  Note  that  both  Che  Hg-rich  and  the  Te-rich  legs  agree  with  experiments  and 
both  of  them  confirm  Che  same  value  ASq  *  9.2  e.u.  The  fact  that  9.2  e.u.  is  also  the 
independently  measured  experimental  value  of  ASq  in  Reference  10  supports  the  validity  of 
the  association  model  and  Che  consistancy  of  the  theory,  including  Vieland's  method  of 
estimating  the  free  energy  of  the  solid  state. 

Figure  A  it  for  the  Cd-Te  binary  liquidus.  The  points  are  experiments  in  Reference  7 
and  the  curve  is  the  present  theory.  The  melting  ceiqMrature  is  fixed  at  Che  experimental 
value  Tq  “^1082*C.  Other  parameters  are  ASn  »  S.O  e.u.,  £,,  “  -1000/6  cal.,  £[3/612  “ 

-5.5,  *23/612  “  "3.3,  £[4  *  *24  "  634  ”  -1000*,  and  u  *  -6000*.  Figure  5  corresponds  to 
Fig.  3  and  plots  the  amount  of  Cd  in  Che  Te-rich  leg  of  the  liquidus  on  a  semi-log  scale. 
The  experimental  points  are  due  to  Harman  (>2).  As  was  discussed  above,  the  slope  of  the 
line  in  Fig.  5  determines  the  value  of  ASn,  which  we  found  Co  be  5.0  e.u.  This  value  of 
ASm  is  different  from  8.8  e.u.  reported  in  Reference  10.  Since  our  estimate  is  regarded 
fairly  reliable,  further  independent  experimental  measurements  of  ASq  for  CdTe  are 
desirable. 
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FIG.  5 

Amount  of  Cd  in  the  Te-rich  side 
liquidus. 


Cd-Te  binary  liquidus.  Points 
are  experiments  (7);  curve  is 
present  theory. 
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Although  Che  aaeociecion  model  la  acill  concroveralal  (13),  the  reaaonably  good 
agreement  with  experimenca  aupporta  ita  uaefulneaa  aa  a  working  model. 

The  pair  approximation  of  CVM  can  take  into  account  the  ahort  range  order,  i.e.,  the 
local  atomic  arrangement,  while  the  point  approximation  cannot.  The  former  ia  recomnended 
in  treacing  the  liquid  acace. 


Acknowledgment 

Diacuaaiona  with  John  W.  Cahn  of  National  Bureau  of  Standarda  and  Lloyd  H.  DeVaux  of 
Hughea  Reaearch  Laboracoriea  are  cordially  acknowledged.  Thanka  are  due  to  Kevin  J.  Riley 
of  Santa  Barbara  Reaearch  Center  for  allowing  the  paper  to  uae  hia  reaulca  before 
publication. 

Referencea 

1.  G.8.  Stringfellow  and  P.E.  Greene,  J.  Phys.  Chem.  Solida  30,  1779  (1969). 

2.  R.  Kikuchi,  Phyaica  103B,  41  (1981). 

3.  R.  Kikuchi,  J.  Chem.  Phya.  60,  1070  (1971). 

4.  R.  Kikuchi,  J.  da  Phya.  (Paria)  38,  C7-307  (1977). 

5.  E.A.  Guggenheim,  Hixturaa  (Oxford  Univeraity  Preaa,  1952). 

6.  L.J.  Vie  land,  Acta  Met.  U_,  137  (1963). 

7.  J.  Staininger,  A.J.  Strauaa,  and  R.F.  Brebrick,  J.  Eleccrochera.  Soc.  117,  1305 
(1970). 

8.  D.A.  Goodman,  J.W.  Cahn,  and  L.H.  Bennett,  to  be  publiahed. 

9.  T.C.  Harman,  "Phyaica  and  Chemietry  of  II-VI  Compounda"  (N.  Aven  and  J.S.  Prener, 
eda.)  p.  767,  North  Holland  Publiahere,  Amacerdam  (1967). 

10.  A.  Laugier,  Rev.  Phya.  Appl.  £,  259  (1973). 

11.  K.J.  Riley,  unpubliahad  (1981). 

12.  T.C.  Harman,  J.  Electronic  Materia la  £,  945  (1980). 

13.  M.  Hillert  and  L.-t.  Staffanaaon,  Met.  Trana.  6B,  37  (1975). 


APPENDIX  C 


THEORETICAL  CALCULATION  OP  Hg-Cd -Te 


LIQUIDUS-SOLIDUS  PHASE  DIAGRAM 


55 


Theoretical  calculation  of  Hg-Cd-Te  liquidus-solidus  phase  diagram  m) 
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The  liquidus-solidus  phase  diagram  of  Hg-Cd-Te  is  calculated  using  the  pair  approximation  of 
the  cluster  variation  method.  The  work  is  an  extension  of  previous  III-V  work  and  includes 
association  as  an  additional  feature.  The  liquid  phase  is  approximated  by  a  pseudo-lattice 
structure  and  includes  molecular  species  as  well  as  atomic  species.  Vacancies  are  also  considered 
in  order  to  take  the  pressure  effect  into  account.  Unlike  the  regular  solution  model,  the  pair 
method  can  take  into  account  the  short-range  order  in  the  liquid  phase.  Adjustable  parameters 
(mostly  energies)  are  first  determined  to  fit  Hg-Te  and  Cd-Te  binary  experiments.  The 
association  model  can  explain  the  asymmetry  and  the  sharp  peak  at  50%  composition.  These 
parameters  and  several  additional  ones  constitute  the  one  set  of  parameters  (independent  of 
temperature  and  composition)  which  are  used  to  calculate  the  entire  ternary  liquidus-solidus 
diagram.  Each  parameter  controls  certain  features  of  the  diagram.  Although  the  work  is  still  in 
progress,  reasonably  good  agreement  with  known  experiments  has  been  achieved  for  both  the  Hg 
comer  and  the  Te  comer.  The  chemical  potential  diagram  is  discussed. 

PACS  numbers:  64.70.Dw,  68.45.  —  v,  81.30.Bx 


Ternary  phase  diagrams  of  III-V  semiconductors  were  pre¬ 
viously  calculated1  using  the  pair  approximation  of  the  clus¬ 
ter  variation  method  (CVM).2  The  work  in  Ref.  1  gives  the 
same  results  as  those  obtained  by  Stringfellow  and  Greene,3 
who  used  the  quasichemical  approximation  of  Guggen¬ 
heim.4  The  present  report  is  an  extension  of  Ref.  1,  and  in¬ 
cludes  association3  and  vacancies  as  additional  features.  Dis¬ 
cussion  of  associated  binary  liquids  can  also  be  found  in  Ref. 
6.  Since  the  theory  is  closely  related  to  Ref.  1,  the  reader  is 
advised  to  consult  Refs.  1  and  6  for  additional  and  pertinent 
background  information.  Compared  with  the  previously 
used  regular  solution  model,5  the  advantage  of  the  pair 
method  is  that  the  short-range  order  in  the  liquid  phase  can 
be  taken  into  account. 

In  the  liquid  phase,  x,  indicates  the  probability  of  finding 
an  ith  species  on  a  lattice  point.  Since  we  will  discuss  both 
binary  and  ternary  systems,  we  will  make  /  =  1,  2,...,  6  de¬ 
note  Hg,  Cd,  Te,  HgTe,  CdTe,  and  a  vacancy  throughout  the 
paper.  For  a  Hg-Te  binary  liquid,  we  use  i  =  1,  3,  4,  and  6 
only,  and  for  a  Cd-Te  binary  we  use  /  =  2,  3,  5,  and  6  only. 
Along  with  the  x,  ’s  we  use  the  probabilities  yv,  for  finding  i 
and  j  species  on  nearest-neighboring  lattice  sites.  These  y/s 
are  the  basic  variables  in  the  pair  approximation  of  the 
CVM.  The  two  sets  of  variables  are  connected  by  the  geo¬ 
metric  reduction  equation 

*<  *  (1) 

The  equilibrium  state  of  the  liquid  phase  is  determined  by 
writing  the  Helmholtz  free  energy  F  in  terms  of  y/s,  and 
then  minimizing  the  grand  potential  G  • sF  —  2,  fs,Nt  with 
respect  to  yv’ s.  (We  write  G  for  the  grand  potential,  not  the 
Gibbs  free  energy.)  The  quantity  p,  is  the  chemical  potential 
(with  respect  to  an  appropriately  chosen  reference  state),  and 
N,  is  the  total  number  of  the  fth  species  in  the  system.  When 
G  is  minimized,  we  obtain  a  set  of  equations 


yv  =  (x(xJ)|I-"/IXexp  -0ev  +  -£-(#*,  + 

(2) 

where z  is  the  coordination  number,  and)?  =  (A7')—  1  and  etJ 
are  the  energy  parameters  for  an  i-j  pair.  We  define  etj  in 
such  a  way  that  eu  —  0  for  the  same  subscripts.  One  advan¬ 
tage  of  the  CVM  formulation  is  that  Eqs.  ( 1 )  and  (2)  hold  for  a 
ternary  as  well  as  a  binary  system,  and  also  for  associated  as 
well  as  nonassociated  liquids  by  appropriately  identifying 
the  subscripts.  Derivation  and  discussions  of  Eq.  (2)  are  in 
Ref.  1 .  In  Eq.  (2),  A.  is  the  normalization  constant  and  has  the 
physical  meaning  of  A  =  —  pv,  where p  is  the  pressure  and  v 
is  the  volume  per  lattice  point. 

The  mathematical  procedure  of  the  problem  is  to  solve 
y/s  from  the  simultaneous  nonlinear  Eq.  (2),  together  with 
Eq.  ( 1 )  for  fixed  values  of  ,  p„  0,  and  z.  They  are  solved  by 
using  the  natural  iteration  method  (NIM)7  without  further 
analytical  transformations.  In  the  NIM,  the  first  input  is  a 
set  of  guess  values  of  x,.  We  use  these  values  on  the  right- 
hand  side  of  Eq.  (2)  to  calculate  yv  as  the  output.  These  y/s 
are  used  in  Eq.  (1)  to  obtain  the  next  input  set,  x<.  It  was 
proved7  that  at  each  iteration  cycle  the  G  value  decreases, 
and  hence  the  iteration  procedure  always  converges. 

In  deriving  the  liquidus  curve  theoretically,  we  need  the 
relation  between  the  free  energies  of  the  liquid  and  solid 
phases.  The  chemical  potential,  for  solid  HgTe  (or 
CdTe)  is  derived  from^,,)  in  the  supercooled  liquid  state  of 
the  50%  composition  by  using  Vieland’s  formula,*  ex¬ 
pressed  as 

fi"'=f*'n-(Tm-T)ASm,  (3) 

where  T„  is  the  melting  point  and  ASm  is  the  entropy  of 
melting  per  molecule. 

For  the  Hg-Te  binary  system,  the  important  adjustable 
parameters  are  el3,  el4,  e)4,  ASm  and  Ut,  where  U,  controls 
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PlC.  1.  Cd-Te  binary  liquidus.  The  curve  is  the  present  theory,  and  the 
points  are  experiments  shown  in  Ref.  9. 


ferent  from  the  III-V  cases,  are  the  asymmetry  and  the  sharp 
peak  at  the  50%  composition.  The  theory  can  derive  asym¬ 
metry  when  The  atomic  pair  interaction  e,3  is  at¬ 

tractive  (i.e.,  negative).  When  f14  and  are  repulsive  (i.e., 
positive),  the  liquidus  shows  a  peak  at  the  50%  composi¬ 
tion.6  When  —  Ux  is  large,  the  association  property  in¬ 
creases  and  the  peak  becomes  sharper.  Then  when  —  f  |«/f  13 
is  larger,  the  left  shoulder  of  the  Hg-Te  liquidus  goes  up,  and 
when  —  e13  is  smaller,  the  entire  liquidus  curve  becomes 
lower  (except  that  Tm  at  the  50%  composition  is  anchored). 
The  value  of  ASm  controls  the  shape  of  the  leg  parts  of  the 
liquidus  curves  in  Figs.  1  and  2. 

The  values  of  the  parameters  we  used  in  calculating  the 
curves  in  Figs.  1  and  2  are:  e,,  =  —  200/6  cal, 
€  1 4/ f  13  a  —  8,0,  ^34/^13  =  “0.5,  U  |  =  —  3000  , 
^23  s  1000/6  cal,  625/^23  s  5.5,  £35/^23  =  —  3.3, 
U2  =  —  6000*,  4£m(HgTe)  =  9.2  e.u.,  Tm(HgT)  =  673  *C, 
ASm  (CdTe)  *  5.0  e.u,  Tm(CdTe)  =  1082  *C,  e*  =  -  1000*, 

i  *  1,  2 . 5,  and  =  0.  The  energy  units  are  left  as  they 

were  used  in  computations.  The  value  of  z  is  chosen  as  6  in 
the  numerical  work;  however,  z  =  4  or  8  does  not  change  the 
results  very  much. 

In  calculating  binary  liquidus  curves,  several  cases  were 
noticed  in  which  different  combinations  of  some  parameters 
led  to  practically  the  same  shape.  Therefore,  the  above  val¬ 
ues  of  the  parameters  are  not  altogether  unique. 

However,  the  value  of  ASm  (HgTe)  is  believed  to  be  accu¬ 
rate  because  9.2  e.u.  leads  to  a  good  agreement  in  both  Hg- 


the  amount  of  association  and  connects  fi,’  s  as 
/*■  +M3-/*4=  Uf  (4) 

Details  of  this  equation  are  found  in  Ref.  6.  Figures  1  and  2 
compare  the  theory  (solid  curves)  of  binary  liquidus  with 
experiments.  Two  features  of  liquidus  curves,  which  are  dif- 


Fic.  2.  Hg-Te  binary  liquidus,  showing  the  present  theory  (curve)  and  ex¬ 
periment*  in  Ref.  10  (points). 


FRACTION  OF  Te 


Fig.  3.  The  Hg  comer  of  the  ternary  Hg-Cd-Te  liquidus.  The  solid  curve  is 
the  present  theory,  and  the  broken  curves  are  Riley’s  experiments.  Ref.  12. 
The  numbers  0.2  and  0.4  are  the  values  of  x  in  the  solid  phase 
(HgTe),  .(CdTe),. 
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FRACTION  OF  Cd 

Fig.  4.  The  ternary  Hg-Cd-Te  liquidus  near  the  Te  comer.  The  solid  curves 
are  the  present  theory,  and  the  broken  curves  are  the  experiments  of  Har¬ 
man  in  Ref.  1 3.  The  numbers  0.2  and  0.4  have  the  same  meaning  as  in  Fig.  3. 


rich  and  Tc-rich  leg  parts  and  also  is  the  value  reported  by 
Laugier." 

For  the  ternary  system  (as  was  done  in  Ref.  1)  we  need  f, , 
the  interaction  of  a  Hg-Cd  pair  on  one  of  two  fee  sublattices 
which  constitute  the  zincblende  structure.  Besides  e,,  ad¬ 
justable  parameters  are  fl2  for  the  Hg-Cd  pair,  f13  for  the 
Hg-CdTe  pair,  e4J  for  the  HgTe-CdTe  pair,  and  e24  for  the 
Cd-HgTe  pair. 

In  comparing  with  other  experiments,  we  first  examined 
the  Hg  corner.  The  value  of  e,  adjusts  the  general  shape  of 
the  curves,  and  the  best  choice  we  could  make  was  e,  =  0. 
The  values  of  e12  and  e„  were  adjusted  to  bring  the  Cd  con¬ 
centration  to  the  right  range.  The  values  we  chose  to  plot  in 
Fig.  3  were  e,2  =  500*  and  e„  =  900T.  For  these  energy  val¬ 
ues,  Cd  species  exist  more  or  less  equally  in  the  atomic  Cd 


Cd/Hg 


Fig.  3.  The  Hg-Cd-Te  liquidus  and  solidus  in  the  Te-rich  side.  This  is  a 
replot  of  Fig.  4.  The  experimental  points  are  those  of  Harman,  in  Ref.  1 3. 
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Fig.  6.  The  Gibbs  diagram  of  the  Hg-Cd-Te  liquidus  and  solidus.  The 
upper  Te-rich  side  is  for  430  *C,  and  the  lower  part  is  for  400  "C.  Note  the 
phase  separation  in  the  lower  part. 


form,  and  the  molecular  CdTe  form. 

After  these  values  are  determined  we  move  to  the  Te-cor- 
ner  in  Fig.  4.  The  shape  of  the  curve  is  adjusted  by  In  the 
region  of  Fig.  4,  HgTe  predominates  over  CdTe;  thus,  when 
e«s  is  attractive,  CdTe  increases  (the  contribution  of  atomic 
Cd  is  negligible).  We  chose  the  value  f43  =  —  120*  to  draw 
Fig.  4.  Figure  3  replots  the  data  in  Fig.  4  in  a  different  way. 
The  circles  are  Hannan’s  experiments13  for  temperatures 
ranging  between  45 0*  and  530  *C.  Although  the  experimen¬ 
tal  points  of  different  temperatures  lie  on  one  curve,  and  the 
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Fig.  7.  The  chemical  potential  diagram  corresponding  to  Fig.  6. 
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600  *C  theoretical  curve  fits  the  experiments  well,  the  theo¬ 
retical  curves  for  other  temperatures  are  shifted. 

It  should  be  pointed  out  that  the  values  of  e, ,  el2,  and 
which  have  been  determined  in  Fig.  3,  have  a  strong  effect  on 
Fig.  4.  The  fact  that  the  same  set  of  values  of  e,,  e,2,  and  *,5 
leads  to  a  satisfactory  agreement  with  experiments  in  both 
Figs.  3  and  4  can  be  counted  as  a  merit  of  the  theory.  The 
value  of  e4S  determined  from  Fig.  4  has  practically  no  effect 
on  Fig.  3  because  the  amount  of  HgTe  in  the  Fig.  3  region  is 
much  smaller  than  the  amount  of  Hg.  Since  no  detailed  ex¬ 
perimental  data  were  available  for  the  Cd-rich  region,  we  left 
e24  undetermined  for  the  Cd-HgTe  interaction  and  used 
*24  =  0  in  computation. 

Figure  6  is  the  Gibbs  diagram  of  the  liquidus-solidus. 
Note  especially  the  tie  line  which  goes  all  the  way  across 
from  the  Hg  corner  to  the  point  near  CdTe,  in  agreement 
with  the  report  of  Ref.  14.  Since  our  value  of  *  ,2  =  500* 
shows  the  repulsion  between  Hg  and  Cd,  there  is  a  phase 
separation  in  the  liquid  phase.  In  Fig.  6  the  state  of 
Hg  ==  0.976,  Cd  si  0.024  coexists  with  the  state  of 
Hg  =  0.027,  Cd  =  0.968.  (These  values  are  not  reliable  be¬ 
cause  *24  =  0  is  tentative.} 

Figure  7  shows  the  chemical  potential  diagram  as  it  corre¬ 
sponds  to  Fig.  6.  The  chemical  potentials  for  Hg,  Cd,  and  Te 
are  adjusted  (linearly  shifted)  so  that  the  sum  of  the  three  is 
zero.  In  that  case,  and  js3  can  be  plotted1  in  the  star 
diagram  of  Fig.  7.  The  upper  curve  is  for  the  Te  comer  and 
the  lower  curve  is  for  the  Hg-Cd  side.  Each  point  on  the  two 
curves  represents  the  coexisting  states  of  the  liquid  and  the 
solid.  A  number  next  to  a  mark  on  the  curves  is  the  composi¬ 
tion  x  in  the  solid  phase,  (HgTe),  _  „  (CdTe)* .  The  point  to  be 
noted  is  that  a  tie  line  in  Fig.  6  is  perpendicular1  to  the  n 
curve  at  the  corresponding  x  point.  The  phase  separation 
near  the  Cd-Hg  edge  in  Fig.  5  is  represented  by  the  discon¬ 
tinuous  change  of  slope  in  the  bottom  part  of  Fig.  7. 

The  pair  approximation  has  been  known  for  many  years, 
and  is  known  to  give  more  reliable  results  than  the  regular 


solution  method.  The  difficulty,  however,  in  applying  the 
pair  treatment  to  a  system  of  many  components  was  that  a 
large  number  of  variables  had  to  be  calculated  by  solving 
nonlinear  equations.  This  difficulty  was  solved  by  the  NIM.7 
Different  from  the  Newton-Raphson  method,  the  NIM 
does  not  need  matrix  inversion,  and  it  guarantees  positivity 
of  all  probability  variables  being  calculated.  The  large  num¬ 
ber  of  variables  does  not  cause  any  problem. 

Although  the  association  model  used  in  this  paper  is  still 
controversial,19  the  good  agreement  between  experiments 
and  the  present  theory,  which  uses  one  set  of  parameters, 
suggests  the  usefulness  of  the  model  as  a  working  hypothesis. 
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continued  interest  in  this  work,  and  also  to  Kevin  Riley  of 
Santa  Barbara  Research  Center  for  allowing  the  paper  to  use 
his  experimental  results  before  publication. 
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APPENDIX  D 


THEORY  OF  SOLID  HgCdTe  WITH  DEFECTS 


1 


INTRODUCTION 


In  the  liquidus  calculations  of  III-V  and  II-VI 
semiconductors  in  References  1/  2  and  3  by  the  author,  and  also 
in  those  by  Brebrick  et  al.4,  the  solid  phase  was  treated  as 
stoichimetric  without  lattice  defects.  In  actuality,  lattice 
defects  are  important  since  the  electrical  properties  of  the 
semiconductor  are  affected  by  the  defects.  In  this  appendix  we 
formulate  the  Hg-Cd-Te  solid  phase,  including  lattice  defects. 

The  solid  Hg-Cd-Te  lattice  is  shown  in  Figure  1 .  It  is 
made  of  two  interpenetrating  fee  sublattices.  On  the  white  fee 
sublattice,  Hg  and  Cd  atoms  are  distributed,  while  the  black  fee 
sublattice  is  occupied  by  Te  atoms.  For  this  examination  we 
allow  vacancies  in  both  sublattices,  and  also  we  include 


antiatoms,  which  means  Te  atoms  on  the  HgCd  sublattice  and  Hg  or 
Cd  atoms  on  the  Te  sublattice. 


In  calculating  the  free  energy,  we  use  the  pair  approxi¬ 
mation  of  the  cluster  variation  method  (CVM) 5.  Although  the  two 
sublattices  are  both  of  the  fee  structure,  each  fee  lattice  is 
disordered  by  itself,  and  hence  there  is  no  frustration  effect 
which  is  characteristic  in  the  ordered  fee  structure. 

Therefore,  the  pair  approximation  is  expected  to  be  a  reliable 
method. 

In  contrast  to  the  ordinary  pair  approximation,  we  use  two 
kinds  of  pairs,  intrasublattice  and  intersublattice;  thus  the 
entropy  expression  of  the  present  formulation  has  not  been  used 
before. 


2 


VARIABLES 


We  use  the  designation  I  and  II  for  the  HgCd  sublattice  and 
the  Te  sublattice.  On  each  lattice  point  we  have  either  a 
Hg,  Cd,  or  Te  atom  or  a  vacancy;  these  four  species  are 
designated  by  i  *  1,2,3  and  4. 


th 


The  probability  of  finding  an  i  species  on  the  sublattice 


I  is  written  as  x  .  .  This  is  a  "point"  variable.  The 


probability  of  finding  an  i-j  nearest-neighbor  pair  in  the  I 
sublattice  is  written  as  yT  .  .;  this  is  the  variable  for  an 
intrasublattice  pair.  The  corresponding  variables,  xTT  .  ancj 

11  f  1 

yTT  . .  are  similarly  defined  for  the  sublattice  II.-  Another 

iA'13  .  .  .  th 

pair  variable,  z..,  is  for  the  probability  of  finding  an  l 


*3  th 

species  on  the  sublattice  I  and  a  j  species  on  a  nearest- 


neighboring  II  sublattice  (note  the  order  of  the  two 


subscripts);  z^j  is  an  intersublattice  variable.  Note  that  z^j 


and  z..  are  different,  while  yT  .  .  *  yT  ... 

J1  A  »  A3  A  '  3  1 

These  variables  are  related  by  the  following  relations: 


X  =  ^ 

XI,i  T 

y  ..  ,E  z  , 

yI,ij  j  zi j  ' 

(2.1a) 

II, i  j 

y  a  E  2 

yH,ij  j  2ji  ’ 

(2.1b) 

The  normalization  relations  are 


1  =  E  y  .  .  a  E.  y  ,  , 
i , j  yI,i3  C3  yII,i] 


-E  ... 

i,j  A3 


(2.2) 


In  subsequent  sections  we  minimize  the  grand  potential 
under  constraints.  The  constraints  written  in  Equation  (2.1) 
can  be  treated  using  Lagrange  multipliers,  o  .  and  <*  .  as 
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GRAND  POTENTIAL 


The  entropy  expression  can  be  derived  by  using,  for 
example.  Barker's  procedure6: 


-6 


(3.1  ) 


where  k  is  the  Boltzman  constant,  N  is  the  number  of  lattice 
points  in  one  of  the  fee  sublattices,  andj£(x)  is  defined  as 

<=C(x)  =  x*nx  -  x  .  (3.2) 


In  writing  the  energy  expression,  we  need  two  kinds  of 
energy  parameters:  for  an  intrasublattice  pair  and  for 

an  intersublattice  pair.  Using  these  parameters,  we  can  write 
the  total  energy  of  the  system  as 


E/N  =  6 


L 

1»  j 


ij 


[yi,ij +  yn.ij] 


+  4 


u 


13  13 


(3.3) 


We  assume  that  a  vacancy  does  not  contribute  potential 
energy,  so  that 


ei4  "  E4i  *  ui4  =  u4i  =  0  for  1  =  1 r 2 , 3 , 4  (3.4) 

The  grand  potential,  8,  is  defined  as 

8  2  E  -  TS  -  5  UiNi  ,  (3.5) 


where  is  the  chemical  potential  and  N^  is  the  total  number 
(both  of  the  i***1  species).  For  the  vacancy  i  =  4,  we 
define  Pi,  *  0,  since  we  do  not  need  the  chemical  potential  for  a 
vacancy.  Using  the  variables  xT  .  and  xTT  . ,  we  can  write 

X  » 1  11 r  1 
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Ni  *  N  (xlti  +  xIIfi) 


(3.6) 


When  we  combine  the  two  terms,  E  -ZujNi  in  Equation 
(3.5), it  is  convenient  to  define 


/e\  .  s  e .  .  -  (e .  .  +  e..  )/2. 
13  13  v  ii  33  " 

u\j  2  u..  -  (uu  +  u..)/2. 


^i  5  wi  -  2uu  "  6  eii  v 


(3.7a) 


(3.7b) 


Then  we  can  write 


5  6  fl  6ij 


[E  -  ?  WiNiJ/N 

[yT,ij +  yn,ij] +  4  h  "u  zij 


i  “i  (xI,i  +  xII,i^' 


(3.8) 


where  each  summation  goes  over  i  =  1,2, 3, 4  and  j  =  1,2, 3, 4. 
From  Equation  (3.7),  the  parameters  with  carets  have  the 
properties: 

^ii  =  “ii  ’  0  and  V  0  1 


(3.9a) 


*i4  =  *4i  =  '  eii/2 


ui4  ■  u4i  -  -  u4,./2 


(3.9b) 


It  is  convenient  to  use  these  energy  parameters  with  carets 
because  in  the  limit  of  no  vacanies,  the  individual  and  u^ 
do  not  contribute  to  the  state  of  the  system. 
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Combining  the  energy  and  chemical  potential  terms  in 
Equation  (3.8)  and  the  entropy  expression  in  Equation  (3.1),  the 
grand  potential,  0,  in  Equation  (3.5),  together  with  the 
constraint  terms,  is  written  explicitly  as  follows: 


• 5  Bti/N  - 66  Ki: +  »n.u] +  46  ft « 

'  Ki + 


/N  /N 
”  .  .  Z  •  • 

13  13 


-  15 


6  [£j‘e(yi,ij)  +^Qyn,ij)] +  4  -  ,4 

-  £ 

4  6  ift  K,j  +  aI,j^  YI,ij  +  6^“ll,i  +  aII,j)  YII,ij 


+  6XZ  (,  -tE  ,  )  ♦  nz  (1  -  ft  yi(..}  +  MlI  u  -  ,r.  yII(iJ 


—  12  fa  +  a  .  )  7 .  .  . 

i,j  ^  i,i  +  ii, y  zij  ' 


(3.10) 


where  the  1  terms  are  for  the  normalizations  and  the  a  terms  are 
the  constraints  derived  in  Equation  (2.3). 
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MINIMIZATION  OF  THE  GRAND  POTENTIAL 


The  equilibrium  state  of  the  system  is  obtained  as  a 
minimum  of  the  grand  potential,  ft,  or  *  in  Equation  (3.10).  In 
minimizing  it,  we  regard  yT  ••,  yTT  and  zij.  as  independent. 

IfiJ  11,  XJ  1J 

We  treat  Xj  ^  and  Xjj  ^  as  averages  of  y's  and  z's;  the  weights 
are  arbitrary,  but  we  use  the  following  expressions: 


(4.1 


Dif ferentrations  of  *  in  Equation  (3.10)  with  respect  to 
yi:i  and  z ^  lead  to  the  following  equations: 


*>  yi,ij  ■  SV6  *  tn  ^i,ij 

ln  yu,u  -  Sin/6  +tn^n,ij 


(4.2 


/\ 


An  z t j  =  0^z/4  +  An  Zij  , 


where  we  separated  out  the  normalization  factors,  x's,  and  the 
quantities  with  carets  are  defined  as 


tn  *I,lj  ‘  it1"  (xx,i  xi,j)  +  6  [‘  +  T6<“i  +  “j>]- 

<n^n,ij  *  if*n(xii,ixii,j)+  s[-^ij  *  +  /“j)]-  0ii,r“xi, 

‘Aj  =d*n(xIfixIIfj)  +  Sp3ij  +  yJCu'i  +  '»j)]  +  3  K,1  *  «u,j 


(4.3 


The  equilibrium  state  is  solved  from  Equation  (4.2)  and 
(4.3),  and  the  constraint  relations  in  Equations  (2.1)  and 


i 


When  these  equations  are  solved,  *  in  Equation  (3.10)  can 
be  simplified  as 


*  S  «  -  .S  y  -  _  2  y  ,  — - 

1/3  yI/iD  3yifij  i/3  yIl/i3  3y 


II, ij  i:>  3zij 


(4.4) 


>  ; 


=  e  (xz  +  Xj  +  xXI)  =  en/N 


On  the  other  hand,  since  the  term  in  Equation  (3.5)  is  the 

Gibbs  free  energy,  we  can  identify  0  as  the  difference  between 
the  Helmholtz  and  Gibbs  free  energies,  so  that 


Q  =  -pv 


(4.5) 


where  V  is  the  total  volume  of  the  system.  From  the  last  two 
equations,  we  can  derive  the  pressure  expression  as 


P  -  -  (*z  +  +  *ir)  N/v 


(4.6) 


Since  we  will  need  it  later,  let  us  derive  the  chemical 
potential  expressions  from  Equation  (4.3)  by  eliminating  o's. 
Making  i  ■  j,  we  obtain 

“1  *  CT[-  J2  tn  (*i,i  *11, i)  +  3  tn  yii.a)  +  2  «*>  2ii] 


”  2  ^z  +  XI  +  *11  ^ 


(4.7) 


This  holds  for  I  ■  1,2, 3, or  4.  The  X  terms  can  be  rewritten  from 
the  i  =  4  case  as 

|(»,+VXII)  ’  kT[-  I2‘n  l*!,!*!!,*'  +  3tn  ^1, 44*11, 44 ]  +2  *"  *44} 


(4.8) 
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5.  ITERATIVE  SOLUTION  OF  THE  EQUATIONS 

The  set  of  equations  in  the  previous  section  is  solved 
iteratively.  The  method  is  called  the  Natural  Iteration 
Method.7  The  flow  chart  of  the  iteration  procedure  is  shown  in 
Figure  2.  In  using  this  scheme  we  fix  the  values  of  the 

.  A  A  .  ,  A 

energies,  e..  the  chemical  potentials,  ui,  and  the 

i  j  /  *  J 

temperature,  3.  The  three  chemical  potential  values  control  the 

equilibrium  composition  and  the  pressure. 

In  the  flow  chart  of  Figure  2,  the  second  step  is  called 

the  Minor  Iteration.  It  solves  the  Lagrange  multipliers  (a's) 

using  the  following  procedure.  The  variable  a  .  is  contained 

J-/1 

in  both  the  yj  terms  and  z  terms  of  Equation  (2.1a),  as  is  seen 
in  Equation  (4.3).  Using  (4.2)  in  (2.1a),  we  write  the  latter 
as 


X _/6  -aT  . 

1  e  1  ? 


°I,i  /\ 
e  ’  v_ 


r*  \  i  • 


(5.1) 

Since  and  are  individually  normalized  to  unity,  one  of 

the  four  equations,  i  =  1,2, 3, 4,  in  Equation  (2.1a)  is 
redundant.  This  means  that  we  need  only  three  otj's.  Or  we 
may  choose,  for  example. 


“l,3  -  0 


a 


11,3 


(5.2) 


Then  we  can  eliminate  the  normalization  factors  X's  in  Equation 
(5.1)  and  write 


MINOR  ITERATION  SOLVES 
aU  &  aD.  i 


When  we  use  this  expression  as  the  iterative  equation  for  ^ , 
we  use  for  the  Oj  ^  written  inside  the  summation  signs, 

and  the  outside  <*  .  is  written  as  a (output)  Thus,  we  write 


(output) 

I/i 


(input) 

I,i 


% 


AT 


T  yi,3j 


(5.4a) 


In  this  equation,  we  use  a  damping  factor,  Y.  We  assign  to  Y  a 
certain  constant  value  between  0  and  1  in  the  iteration 
calculation.  When  Y  *  0,  the  output  a  is  equal  to  the  input 
a  and  the  iteration  does  not  proceed.  When  Y  =  1 ,  the 
0(output),  Which  is  the  next  input  value,  may  sometimes 

if  1 

overshoot.  Usually,  Y  =  0.5  leads  to  stable  convergence  of  the 
iteration. 

For  <*TT  .,  the  equation  corresponding  to  Equation  (5.4a)  is 

II  #  X 


(output) 
II, j 


(input) 
II, j 


JL/-UL1.U 
i  yH,i3 


(5.4b) 


The  chemical  potential  controls  the  composition  as  follows. 


We  first  combine 


x_  .  and  xTT  .  and  define 

i  ,  1  XX  ,  X 


X 


i 


+ 


(5.5) 


Using  xi's,  we  define 


n2  s  xx  -  x 2 


n3  2  X1  +  x2  "  X3 


(5.6) 


Then,  h 2  signifies  the  Hg-Cd  composition,  and  n3  represents  the 
deviation  from  the  stoichiometry.  Using  n's,  we  can  write  the 
chemical  potential  terms  in  the  grand  potential  as 
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where 


or 


"u  S  ^  -^2 

5  ^1  +  *2  ~  'V  2 

\  +  ^2  )/ 4  +  ^3/2 


'y'  =  /y'  +  ^  +  /y' 

1  V  w  u 


B  <S  .  -<S  /\ 

2  *v  w  *u 


^3  “ 


(5.8a) 


(5.8b) 


Prom  Equation  (5.7),  we  see  that  the  value  of  controls 
1  -  x4;  as  uv  increases,  1  -  x4  increases  and  the  vacancy 
concentration  decreases.  The  value  lof  controls  the  Hg-Cd 
composition;  as  /fru  increses  xj  -  x2  increases.  As  the  value 
increases,  the  non-stoichiometry  Xl  +  x2  -  x3  increases. 

When  we  treat  a  ternary  system,  we  expect  it  to  reduce  to 
the  binary  case  when  one  of  the  ternary  components  reduces  to 
zero.  For  example,  when  x2  reduces  to  zero,  /mv2  goes  to  -  *, 
whereas  'v'j  and  ^2  remain  finite.  We  can  establish  the 
correspondence  between  the  ternary  and  the  binary  case  by  this 
procedure. 
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SUMMARY 


The  Hg-Cd-Te  crystalline  state  is  formulated  using  the  pair 
approximation  of  the  CVM,  including  vacancies  and  antiatoms,  but 
without  interstitial  atoms.  The  crystal  is  made  of  two  fee 
sublattices  I  and  II.  On  each  lattice  point,  one  of  the  four 
species  is  found:  Hg,  Cd,  and  Te  atoms  and  a  vacancy.  The 
treatment  uses  nearest-neighbor  intersublattice  potentials,  u^j, 
and  nearest-neighbor  intrasublattice  potentials,  e^j.  The  basic 
variables  are  intersublattice  pair  probabilities,  z. . ,  and 
intrasublattice  pair  probabilities,  yT  .  .  and  yTT  ... 

if  ij  Aif  ij 

In  order  to  solve  the  equilibrium  state,  the  grand 
potential,  0  =  E  -  TS  -£i  Ui  Nj.,  is  written  in  terms  of  y's  and 
z's,  and  is  then  minimized  with  respect  to  the  variables.  The 
resulting  equations  are  solved  for  given  values  of  T  and  Mi's  by 
an  iterative  technique. 

This  solid  state  treatment  is  a  part  of  the  work  on  the 
liquidus-solidus  calculation  in  Appendix  E. 
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A.  INTRODUCTION 

The  liquidus  and  solidus  of  the  binary  Hg-Te  and  Cd-Te 
systems  and  the  ternary  Hg-Cd-Te  have  been  calculated  by  the 
author1'2  using  the  pair  approximation  of  the  cluster  varia¬ 
tion  method  (CVM).  In  these  calculations,  the  free  energy,  F, 
of  the  binary  solid  state  was  represented  by  a  point  in  the  F 
vs.  composition  space  because  the  solid  phase  was  assumed  to  be 
stoichiometric  and  free  from  any  lattice  defects.  Since  the 
lattice  defects  are  of  practical  importance,  we  included  vacan¬ 
cies  and  antiatoms  in  the  solid  phase  in  Appendix  D  of  the  pres¬ 
ent  report.  This  appendix  reports  the  procedure  of  combining 
the  new  solid  phase  in  Appendix  D  and  the  liquid  phase  of  Refer¬ 
ence  2  to  obtain  the  liquidus  and  the  solidus. 

In  the  ternary  liquid  treatment  of  Appendix  C,  not  many 
equations  were  written  because  of  the  page  limitation.  Actu¬ 
ally,  the  equations  in  Appendix  B  for  the  binary  liquid  can  be 
ujied  for  the  ternary  liquid  as  well  without  much  change.  Thus, 
w?e  may  refer  to  Appendix  B  when  equations  for  a  ternary  liquid 
are  needed. 

B.  COEXISTENCE  CONDITIONS 

When  two  phases  coexist,  one  of  the  conditions  of  coexist¬ 
ence  is  that  the  chemical  potentials  of  each  component  are  equal 
in  the  two  phases;  i.e., 

MiU)  .  Wl(s)  i  -  1,  2  and  3,  (2.1) 

where  1,  2  and  3  denote  Hg,  Cd  and  Te  species.  It  should  be 
mentioned  here  that  the  previous  work  of  the  Hg-Cd-Te  liquidus 
and  solidus2  could  not  use  Equation  (2.1)  because  the  solid 
phase  was  represented  by  a  singular  structure,  as  was  mentioned 
in  the  previous  section,  and  because  was  indeterminate. 

The  second  condition  of  coexistence  is  that  the  pressure  of 
the  liquid  and  that  of  the  solid  are  the  same.  Since  our  models 
of  the  liquid  and  the  solid  both  contain  vacancies,  we  can  cal¬ 
culate  pressures  of  both  phases  and  equate  them  for  coexistence. 


When  we  compare  chemical  potentials  in  the  liquid  and  solid 
phases,  we  need  to  determine  the  reference  levels  in  the  two 
phases.  We  use  Vieland's  equation3,  as  was  done  in  Reference  2. 
It  is  written  for  binary  cases  as  follows: 

(s,i3)  y^(s,i3)  a  (sci.,i3)  +  y^^scA/iS)  _  (Tj^  —  T)  As^, 

i  =  1  and  2  (2.2) 

where  s  and  sc  A  in  supercript  stand  for  solid  and  supercooled 
liquid,  respectively,  and  i3  indicates  that  this  equation  is  for 
the  i-3  binary  case.  In  the  last  terra,  Tm  is  the  melting 
temperature,  and  ASm  is  the  entropy  of  melting  (per  molecule) 
of  HgTe  (i  =  1)  or  CdTe  (i  *  2)  solid. 


C.  STEPS  OF  LIQUIDUS  -  SOLIDUS  CALCULATION 

The  liquidus-solidus  calculations  are  done  in  the  following 
steps.  The  temperature  T  and  the  pressure  p  are  fixed 
throughout  the  calculation.  In  this  and  subsequent  sections  we 
use  the  names  i  *»  1,  ..,6  for  atomic  Hg,  Cd,  Te,  molecular  HgTe, 
CdTe  and  a  vacancy  in  the  liquid  phase. 

(1)  First,  we  solve  the  Hg-Te  and  Cd-Te  binary  liquid 
states  for  the  50-50  composition.  This  state  is 
unstable  below  the  melting  temperature,  Tm,  and  was 
called  the  supercooled  state  in  Section  2,  with  the 
designation  sc*  in  Eqyation  (2.2).  We  calculate  the 
sum,  m(s,i3)  +  w3(s'l3>,  for  a  binary 

solid  phase  from  the  left-hand  side  of  Vieland's 
relation  (Equation  2.2). 

(2)  Then,  we  solve  the  binary  solid . states,  Hg-Te 
and  Cd-Te,  and  ^(s»i3)  +  ^3(8,13)  are 

obtained.  The  superscript  (s,i3)  indicates  that  this 
is  the  i3  binary  solid  state  (i  -  1  or  2). 

(3)  Next,  we  solve  the  ternary  solid  state.  As  we 
remarked  in  Section  5  of  Appendix  D, 
controls  the  HgsCd  ratio  and  ^w(s)  controls 

the  (Hg  +  Cd):  Te  stoichimetry .  The  third  chemical 
potential  parameter,  controls  the 

number  of  vacancies  and  thus  the  pressure.  By 


selecting  these  three  chemical  potential  parameters, 
we  make  the  solid  phase  take  the  desired  composition 
and  the  prescribed  pressure,  pQ.  For  this  state, 
the  three  values,^!,  in  Equation  (5.8b)  of  Appendix 
D  are  written  as  /Q^(s,123)  to  note  that  it  is 
for  the  ternary  case. 

(4)  In  the  solid  phase,  u^  and  are  related  by  Equation 
(3.7b)  of  Appendix  D: 

/VS)  -  Vi(S)  “  2uu  -  6eu  .  (3.1) 

We  assume  that  and  remain  the  same,  indepen¬ 
dent  of  the  ternary  composition.  Then  we  can  write 

u  .  (s, 1 23)  y  ( s ,  1  23 )  _  ij  ( s ,  i  3 )  j,(s,i3) 

Mi  3  Mi  3 

+  ^.(s,123)  +  ^3(s,123)^ 

_  /^(s,i3)  +  /^(s,i3)J 
i  *  1  and  2  .  (3.2) 


On  the  right-hand  side  the  first  two  terms  are 
calculated  using  the  Vieland  formula  (Equation  2.2), 
and  the  terms  with  carets  are  calculated  in  solving 
the  solid  phases.  Thus,  we  can  evaluate  the  left-hand 
side  quantity. 


(5)  In  the  coexisting  states,  Uj/s  are  equal  in  the 

liquid  and  the  solid  phases,  as  was  stated  in  Equation 

(2.1) .  Thus  we  can  replace  the  u's  on  the  left-hand 
side  of  Equation  (3.2),  by  the  corresponding  u's  for 
the  liquid  phase.  For  the  first  two  terms  of  the 
right-hand  side  of  Equation  (3.2),  we  use  Equation 

(2.2) .  Then  we  can  write. 


+  p3(M  =  Pi(sc*,i3)  +  U3(sc*,i3) 


-  (Tm-T)  as 


+  ^  C  s ,  123)  ^3  ( s ,  123)^ 

-  ^<s' 


m 


■  i3)  +  ^3(s,i3)j 


1  and  2 


(3.3) 
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Since  both  and  yMscA'13)  are  for  the  liquid 
phase,  we  can  use  Equation  (3.b)  of  Appendix  B  and 
place  carets  on  them  to  write 


/fr.  U)  +  ><s(4)  =  ft.  (sc£,i3)  +  /jj,  (sc£,i3) 
+  fo.(s,123)  +  ^3(s,123)J 
_  ^(s,i3)  +^3(s,i3)^ 


-  (Tm-T)  A  Sr 


i  *  1  and  2 


(3.4) 


(6)  We  choose  'fe  j ^ .  Then,  combining  with  Equation  ( 3 . 4 ) , 
we  know'tfj^*  and  ^  •  Using  these  three 
values,  we  can  solve  the  liquid  state.  When  the 
liquid  state  is  solved,  we  examine  the  pressure.  We 
then  vary  to  come  to  the  preassigned  value,  pQ, 

of  the  pressure.  By  calculating  ^  from  'frj ^ '  using 
Equation  (3.6)  and  (3.15)  of  Appendix  B,  we  get 

p3(£)  s  'fjSjf*')  +  u>£23  -  e3  ^  (3.5) 

Note  that  these  are  different  from  the  solid  relation 
in  Equation  (3.1). 


(7) 


Since  we  know  the  value  of  ^3 ( s )  in  Step  (3),  we  can 
calculate  u3Js|  from  Equation  (3.1).  We  then  plot 
p 3 J s '  and  u3{  ^against  ^w(s)  used  in  Step  3.  We  vary 
and  find  the  value  for  which 


(3.6) 


This  point  is  the  coexistence  point.  Note  that  in 
going  from  Equation  (3.2)  to  Equation  (3.3),  we  use 
the  relationship 

ui(S)  +  M3(S)  s  yi<il)  +  y3U)  for  i  *  1  and  2  .  (3.7) 


The  combination  of  Equations  (3.6)  and  (3.7)  are 
equivalent  to  the  coexistence  condition  (Equation  2.1) 


D.  SUMMARY 

We  have  described  the  procedures  of  calculating  the 
liquidus-solidus  of  the  Hg-Cd-Te  system.  For  the  solid  phase, 
the  method  in  Appendix  D  is  used  and  the  lattice  defects 
(vacancies  and  autiatoms)  are  included.  The  liquid  phase  is 
modeled  using  a  pseudo-lattice  structure,  as  was  done  in 
Reference  1  and  2. 

The  detailed  steps,  (1)  through  (7),  of  computing  the 
liquidus  and  solidus  for  the  coexistence  condition  are  described 
in  Section  3. 

The  main  results  are  shown  in  Figures  10,  11  and  12  of  the 
main  text  of  this  report. 


